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1 Executive Summary 

Under the auspices of the Aviation Safety Program (AvSafe) Integrated Intelligent Flight 
Deck (IIFD) project, NASA is investigating remote sensors for detection of hazards to air¬ 
craft. The Forward-Looking Interferometer (FLI) is a new instrument concept for efficiently 
obtaining the measurements required to alert flight crews to potential weather hazards to 
safe flight. To meet the needs of the commercial fleet, such a sensor should address mul¬ 
tiple hazards to warrant the costs of development, certification, installation, training, and 
maintenance. The FLI concept is based on passive radiometric sensing with high-resolution 
Infrared (IR) Fourier Transform Spectrometry (FTS) technologies that have been devel¬ 
oped for satellite remote sensing, in addition to the detection of aerosols and gasses for 
other purposes. The FLI concept is being evaluated for its potential to address multiple 
hazards including clear air turbulence (CAT), volcanic ash, wake vortices, low slant range 
visibility, dry wind shear, and icing, during all phases of flight (takeoff, cruise, and land¬ 
ing). Additional possible applications, not addressed in this study, may include detection 
of runway contamination; observations of temperature and water vapor for use in weather 
models; and detection of pollution or other aerosols. The technology might also be ap¬ 
plied for downward-looking atmospheric profiling of water vapor, temperature, turbulence 
or other applications. 

Over the years, a large suite of techniques and analytical tools has been developed to 
assess the potential measurement performance of downward-looking imaging interferometers 
that were developed for space-based Earth remote sensing. Some of these analytical tools 
have been applied to an airborne hazard detection feasibility study for CAT, volcanic ash, 
and low slant range visibility conducted for NASA by the Georgia Tech Research Institute, 
Hampton University, and the National Center for Atmospheric Research through National 
Institute of Aerospace Task Order Number, NNL06AB99T. In the feasibility study, remote 
sensing techniques developed for space-based sensors were applied to aviation hazard detec¬ 
tion problems in order to assess the feasibility of using FLI technology for airborne sensors. 
Sensitivity studies were performed on clear air turbulence, volcanic ash clouds, and for re¬ 
duced slant range visibility. Characterization studies and literature searches were conducted 
on wake vortices, dry wind shear, and icing to determine if future sensitivity studies were 
warranted. Sensitivity studies were conducted using a limited set of temperature, wind, 
and humidity fields. The studies were performed spectrally with high enough resolution 
to resolve individual spectral lines of molecules such as CO2, so as to support the use of 
interferometric sensors such as the FLI. The simulations were performed for single lines of 
sight from the sensor to/through the hazards. 

The studies were based on a high spectral resolution spherical shell radiative transfer 
model using LBLRTM (Line-By-Line Radiative Transfer Model) for the molecular optical 
depths, and values obtained from the literature for volcanic ash and obscurants such as 
fog and haze. Two different approaches for detecting CAT were considered: detecting 
temperature disturbances ahead of the aircraft with measurements in the 14 - 15 /jm CO2 
spectral band, and detecting water vapor disturbances using wavelengths in the 5-8 micron 
spectral region. For volcanic ash detection, two different algorithms for determining the 
distance to an aerosol cloud were examined. The algorithms are the “CO2 Slicing” technique 
and “Window Channel” technique. 

For the improvement of slant range visibility, the use of Empirical Orthogonal Functions 
(EOF) Regression [also called Principle Components (PC) Regression] was investigated. By 
definition, the EOFs, or PCs, are statistically optimal in the sense that no other set of 
spectral functions can be used to represent the observed radiance spectra with as small a 
number of independent variables. An advantage of using EOF regression is that it enables 
the entire spectrum of information content to be used for the sensing of these phenomena, 
thereby multiplexing the emission signals from all the atmospheric gases of interest as well 



as those emission signals from clouds and aerosols. By using EOF regression to retrieve 
the geophysical characteristics of the otherwise small radiance turbulence signal, the signal- 
to-noise ratio of the FLI detection system may be enhanced by an order of magnitude or 
more. A similar signal-to-noise advantage could be gained for the detection of other weather 
related hazards as well. 

The physical basis for FLI detection of wake vortices, dry wind shear, and icing has 
been investigated through literature research to determine whether sensitivity studies are 
warranted. In all cases, these hazards have phenomenology that should be detectable by an 
FLI system. Wake vortices may be detectable due to entrained water vapor from the aircraft 
engine exhaust, lower than ambient core temperature, and/or other entrained gasses; dry 
wind shear by the cooler temperature of air in a downdraft; and icing by the different optical 
properties of water and ice in the infrared region. 

The initial feasibility studies indicate that the spectral radiance measurements by a FLI 
mounted on an aircraft should be able to sense or, in the case of reduced visibility, mitigate 
aviation weather hazards to safe flight. It was shown that, for the case of enhancing slant 
range visibility, Eigenvector EOF regression is a powerful method for transforming FLI 
spectra into products that can be used by the flight crew or by automated agents to help 
insure safe flight. It is likewise expected that EOF regression can be used to enable rapid 
and accurate predictions of the other aviation safety hazards including turbulence, volcanic 
ash, wake turbulence, and icing. 

Future research needs have been determined to identify the next steps to be taken in 
the development of a FLI-based comprehensive external hazard monitor. Such future ef¬ 
forts should include both ground- and flight-testing of prototype FLI hardware to collect 
data for sensor and phenomenological characterization and validation, in addition to the 
development of means for measuring hazard intensity or estimating the hazard posed to the 
aircraft and the development of a simulation environment. The next modeling steps in the 
FLI applications research should be to develop and apply the EOF regression technique to 
predicting the intensity and distance / warning time for turbulence, volcanic ash, wake tur¬ 
bulence, and icing. The technique should be demonstrated for a wide variety of atmospheric 
situations to demonstrate the robustness of the algorithm. 

CAT and Dry Wind Shear. Sensitivity studies for CAT detection included two tech¬ 
niques: sensing temperature fluctuations using the 650 - 780 cm -1 (wavenumbers 6.5 - 7.8 
m -1 ) CO2 band and sensing water vapor fluctuations in the 1300 - 1600 cm -1 (wavenumbers 
13 - 16 m -1 ) region. The simulations revealed detectable signatures in both the interfero- 
gram and spectral domains for the CO2 temperature technique at flight altitudes. Although 
the maximum detection range appears to be about 5 km (assuming an instrument resolution 
of 0.3 K), which is not as large as early investigators have claimed. But humidity anomalies 
in the one case for which data exists were detectable from 30 km, or with a warning time of 
approximately 2 min, assuming a lower instrument sensitivity of 1.0 K. Wind shear near the 
ground is associated with the downward motion of cooler air from higher altitudes, and the 
anomalously low air temperatures should be easily detectable with an FLI during a landing 
approach. 

Volcanic Ash. The ability of an IR imager to distinguish volcanic ash clouds from 
water and ice clouds was verified in the 800 - 1100 cm -1 region, using an ash optical depth 
spectrum obtained from the literature. The optical properties of volcanic ash produce a 
large signature in the 8-12 cm -1 window region, thereby making ash detection possible 
with the FLI. 


Low Visibility. A simulation of multi-band imaging in low slant-range visibility con¬ 
ditions illustrated the potential of the FLI to provide essentially unobscured IR images of 



runway scenes as viewed from the air during an approach by making use of multi-spectral 
information to remove the effects of atmospheric obscurants (e.g. water vapor and aerosols). 
This capability will give pilots a clear view of runways even when an obscurant such as a 
ground fog is present, and it will provide high quality images to support automated detec¬ 
tion of airborne and ground obstacles or other hazards. 

Wake Vortices. Typical wake vortex properties obtained from a literature search 
included the core temperature as well as the concentrations of exhaust gas constituents 
(including water vapor) that become entrained in the vortex. Vortex core temperatures 
have been measured to be several degrees below ambient temperature and this temperature 
difference should be detectable, which suggests that an FLI will provide real-time imagery 
of the motion and evolution of wake vortices. The entrained gasses, including water vapor, 
provide additional means of detection with an FLI. Detection of wake vortices by FLI tech¬ 
niques is possibly the most challenging application considered due to the relatively small 
size of the phenomena and the geometry of the problem which would have the FLI looking 
down at a scene with ground, runway, and other objects as background during landing, the 
phase of flight that would be of most interest. 

Icing. The optical properties of liquid water and ice are quite different in the spectral 
regions where the FLI will operate. This fact will give the FLI a capability to distinguish ice 
from liquid water on surfaces such as runways. In addition, the FLI may be able to detect 
icing conditions ahead of an aircraft during flight by discriminating supercooled liquid water 
drops from ice crystals. 

Runway Friction. One significant cause of aircraft accidents is reduced runway friction 
caused by rain, frost, and snow. A relatively low spectral resolution imaging FLI spectrom¬ 
eter, operating from an aircraft, could be used to detect this hazardous runway condition. 
Using FLI data, the braking distance can be determined, and used to alert the pilot of a 
hazardous landing condition. 

The sensitivity and characterization studies addressed the phenomenology that supports 
detection and mitigation by the FLI. Techniques for determining the range, and hence warn¬ 
ing time, were demonstrated for several of the hazards, and a table of research instrument 
parameters was developed for investigating all of the hazards discussed above. This work 
supports the feasibility of detecting multiple hazards with an FLI multi-hazard airborne 
sensor, and for producing enhanced IR images in reduced visibility conditions; however, 
further research must be performed to develop a means to estimate the intensities of the 
hazards posed to an aircraft and to develop robust algorithms to relate sensor measurables 
to hazard levels. In addition, validation tests need to be performed with a prototype sys¬ 
tem. Further research should include ground-based measurements, airborne measurements, 
additional sensitivity studies, and simulations. 


2 Introduction 

In the fall of 2005, NASA’s Aeronautics Research Mission Directorate (ARMD) initiated 
a restructuring of its aeronautics program to ensure that it had a strategic plan in place 
that enables the pursuit of long-term, cutting-edge research for the benefit of the broad 
aeronautics community. The three principles guiding this restructuring were: 1) NASA 
will dedicate itself to the mastery and intellectual stewardship of the core competencies of 
aeronautics in all flight regimes; 2) NASA will focus its research in areas appropriate to its 
unique capabilities; and 3) NASA will directly address the fundamental research needs of 
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the Next Generation Air Transportation System (NextGen - formerly known as NGATS) 
while working closely with its agency partners in the Joint Planning and Development Office 
(JPDO). 

Using the above principles, ARMD established four programs: the Fundamental Aero¬ 
nautics Program, the Aviation Safety Program (AvSafe), the Airspace Systems Program, 
and the Aeronautics Test Program. The Fundamental Aeronautics Program conducts 
cutting-edge research that produces innovative concepts, tools, and technologies that en¬ 
able the design of vehicles that fly through any atmosphere at any speed. The Aviation 
Safety Program focuses on developing the cutting-edge tools, methods, and technologies to 
improve the intrinsic safety attributes of current and future aircraft that will operate in the 
NextGen. The Airspace Systems Program directly addresses the Air Traffic Management 
(ATM) research needs of the NextGen initiative as defined by the JPDO. The Aeronautics 
Test Program ensures the strategic availability and accessibility of a critical suite of aero¬ 
nautics test facilities that are deemed necessary to meet aeronautics, agency, and national 
needs. 

Both the restructuring of NASA’s aeronautics research program and the development 
of the National Aeronautics R&D Policy took into consideration many of the same reports 
and studies that have been conducted over the past several years. These reports included 
the “Final Report of the Commission on the Future of the U.S. Aerospace Industry” by the 
Aerospace Commission, 2002, “Securing the Future of U.S. Air Transportation: A System 
in Peril” by the National Research Council (NRC), 2003, and “Wind Tunnel and Propulsion 
Test Facilities: An Assessment of NASA’s Capabilities to Serve National Needs”, RAND 
Corporation, 2003. The use of this supporting information by both parties helped contribute 
to the alignment of NASA’s restructured aeronautics program with the new Aeronautics 
RfeD Policy. In addition, while the recent NRC “Decadal Survey of Civil Aeronautics” 
(2006) was published after the restructuring of NASA’s aeronautics program had begun, 
NASA’s new aeronautics program is very well aligned with the recommendations of that 
study. 

The Aviation Safety Program also conducts research in several critical core aeronautics 
competencies, including materials and structures, sensors, data mining, flight mechanics, 
human factors, verification and validation of complex systems, flight dynamics and control, 
icing physics, and intelligent and adaptive systems. Mastery of these competencies will 
be critical to our ability to develop innovative concepts, tools, methods, and technologies 
that will improve the intrinsic safety attributes of current and future aircraft operating in 
the National Airspace System (NAS). In addition, just as in the Fundamental Aeronautics 
Program, much of the foundational research conducted in the Aviation Safety program will 
ultimately have applicability to spacecraft as well as aircraft. For example, research in core 
areas such as sensors and data mining in the Integrated Vehicle Health Management project 
may lead to self-healing systems that are applicable not only to aircraft but to spacecraft 
as well. Foundational research in materials that leads to advances in our understanding of 
the aging and durability properties of composites and superalloys will have dual aeronautics 
and space applications. Foundational research in adaptive controls as part of the Integrated 
Resilient Aircraft Control project will likely find applications in future space missions where 
vehicles will be required to operate in and adapt to unknown flight environments. Finally, 
foundational research in human factors, sensors, and human-automation interface modeling 
will lead to enhancements in crew-vehicle interface technologies that may have applicability 
to space vehicles as well. 

The Aviation Safety Program worked closely with the FAA and the JPDO in devel¬ 
oping its research plans with a focus on developing cutting-edge technologies to improve 
the intrinsic safety attributes of current and future aircraft that will operate in the global 
NextGen. Furthermore, many of the technologies will have dual-use applications for DOD 
needs. As such, all four projects within the Aviation Safety Program directly support this 



Policy guidance. 

Specifically, the Integrated Intelligent Flight Deck (IIFD) project will pursue flight deck 
related technologies that ensure crew workload and situation awareness are both safely 
optimized and adapted to the future operational environment as envisioned by the NextGen. 
A key component of this research will be investigating methods to automatically monitor, 
measure, and assess the state of crew awareness. Project results should enable system 
designers to eliminate the safety risk of unintended consequences when introducing new and 
advanced systems into an operational environment. 

As stated previously, IIFD develops flight deck technologies that mitigate operator-, 
automation-, and environment-induced hazards for future operational concepts. The scope 
includes: (a) development of crew/vehicle interface technologies that reduce the risk of pilot 
error; (b) development of monitoring technologies to enable detection of unsafe behaviors; 
(c) development of fail-safe methods for changing the operator/automation roles in the 
presence of detected disability states; and (d) development of a comprehensive surveillance 
system design that enables robust detection of external hazards with sufficient time-to-alarm 
for safe maneuvering to avoid the hazards. 

For our purposes, a flight deck and a flight deck system are defined as follows: 

• Flight Deck: A volume of space designed to accommodate at least one human operator 
and the interfaces between the operator and the remainder of the flight deck system. 

• Flight Deck System: A system that includes (1) the entity(s) who have the authority 
and responsibility for directing the flight of an aircraft, (2) all subsystems that directly 
interface to these entity(s), and (3) the interfaces between them. 

Vigilant observation of the surrounding airspace for external hazard detection is essen¬ 
tial for operational capability and safety. There are five classes of external hazards: (1) 
meteorological (e.g. icing conditions, convective weather, wind gusts, turbulence), (2) en¬ 
vironmental (e.g. volcanic ash), (3) geospatial (e.g. terrain, man-made obstacles, foreign 
object debris), (4) traffic, and (5) airspace constraints/restrictions. The severity of hazards 
associated with any of these depends on many factors such as relative proximity and closure 
rate. 

The Forward-Looking Interferometer (FLI) is an instrument concept for efficiently ob¬ 
taining the measurements required to alert flight crews to potential weather hazards to safe 
flight. The FLI concept is based on passive radiometric sensing with an aircraft-mounted 
Michelson Fourier transform spectrometer capable of measuring atmospheric conditions 
ahead of the aircraft for the purpose of detecting hazards to the aircraft or determining 
other information related to flight deck requirements. This instrument concept was initially 
evaluated for detecting Clear Air Turbulence (CAT), and the project was subsequently ex¬ 
panded to include the investigation of the feasibility of detecting other hazards such as 
volcanic ash, wake vortices, dry wind shear, reduced visibility, and icing conditions along 
the flight track of the aircraft. It is potentially applicable to external hazards in classes 
1-4 above. FLI sensor data may be used alone or in combination with other sensor data 
such as airborne weather radar or lidar. The data it provides should complement other 
airborne forward-looking hazard detection, such as radar detection of turbulence or wind 
shear, by providing data in conditions where other sensors have insufficient signal available 
to allow detection. It is also anticipated to be capable of providing additional data useful 
in discrimination and hazard estimation. Additional possible applications, not addressed 
in current investigations, may include detection of runway contamination; observations of 
temperature and water vapor for use in weather models; and detection of pollution or other 
aerosols. The technology could also be applied for downward-looking atmospheric profiling 
of water vapor, temperature, turbulence or other applications. 
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The study seeks to predict instrument performance for hazard detection and measure¬ 
ment via models and analyses, to validate the models with experiments and to identify 
technical challenges that must be addressed for development of a practical sensor based on 
the technology. Subsequent to a successful feasibility study [1], NASA expects to attack 
the technical barriers to instrument development, develop and improve the models and per¬ 
formance predictions, and to identify the next steps to be taken in the development of a 
FLI-based airborne sensor. Such future efforts should include ground and flight-testing of 
prototype FLI hardware to collect data for sensor and phenomenological characterization 
and validation. 

This feasibility study began with a Workshop on the Airborne Radiometric Detection 
of Aviation Hazards. The goals of the workshop were to review the current status of re¬ 
search in this technical area, to assess the resources that are available for further work, and 
to familiarize interested participants with the research program. The workshop included 
participants from government, industry, and academia. The results of the Workshop are 
presented in Appendix A of this report. 

This report includes sensitivity studies and the characterization of various hazards oc¬ 
curring at cruise altitudes or during the takeoff and landing of aircraft. Sensitivity studies 
were conducted on CAT, volcanic ash, and low slant range visibility. Characterization to 
determine if sensitivity studies were warranted included wake vortices, dry wind shear, and 
icing. 

Based on the results of the studies discussed in this report, a table of FLI radiance mea¬ 
surement requirements for a research instrument is presented. In addition, future research 
needs have been determined to identify the next steps to be taken in the development of 
an FLI-based airborne hazard sensor. Such future efforts should include further sensitiv¬ 
ity studies, ground and flight-testing of prototype FLI hardware to collect data for sensor 
and phenomenological characterization and validation, and development of a simulation 
environment. 

Additional information on previous wind shear studies is presented in Appendix B and 
candidate research instruments are described in Appendix C. A report by Larry Cornman 
of the National Center for Atmospheric Sciences (NCAR) is included as Appendix D. The 
NCAR report is divided into two parts: Part A discusses the work performed investigating 
methodologies for measuring turbulence from an airborne IR interferometer; Part B provides 
an overview of the data sets that were delivered to Hampton University in support of their 
simulation activities. 

A feasibility study of using an FLI to detect a hazardous runway surface conditions 
(reduction of runway friction from ice or snow) is the subject of Appendix E, a report by 
William Smith. 


3 Background 

The long history of investigations into airborne IR radiometric detection of aviation hazards 
has provided examples of apparent detections of CAT and wind shear, as well as a data 
set that was interpreted as images of wake vortices. In addition, ground-based radiometric 
detection of volcanic ash clouds has been demonstrated [2], [3], [4]. However, all of the early 
work was done with simple instruments that had low spectral, spatial, and radiometric 
resolution, and most of them employed a single line of site. The detections were generally 
not supported by truth data, and the atmospheric conditions that the radiometers were 
responding to were characterized poorly, if at all. Here we define “truth data” as independent 
measurements of the hazard, e.g. aircraft acceleromteter measurements of turbulence or F- 
factor measurements for wind shear. Because the atmospheric state was not known, the 
performance of the radiometers as hazard sensors could not be checked with numerical 
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simulations. These problems led Gimmestad et al. [3] to recommend that flight tests should 
be conducted with a hyperspectral IR imaging system; that the flights should include truth 
data; and that the measurement effort should be supported by a parallel simulation effort. 
The study reported here is a related effort aimed specifically at assessing the feasibility of 
using an FLI for aviation hazard detection. 


4 Sensitivity Studies 

Sensitivity studies were performed on CAT, volcanic ash, and low slant range visibility. 
CAT sensitivity studies included the development of a high spectral resolution spherical 
shell radiative transfer model, based on Line-By-Line Radiative Transfer Model (LBLRTM) 
[5], [6] molecular optical depths. FLI simulations were conducted using three-dimensional 
temperature and water vapor fields associated with CAT provided by the National Center 
for Atmospheric Research (NCAR) in Boulder, CO. Volcanic ash cloud detection sensitivity 
studies have been conducted using window region spectral radiance emission observations. 
Two different algorithms for determining the distance to an aerosol cloud were implemented. 
The algorithms include the “CO2 Slicing” technique and the “Window Channel” technique. 
Lastly, sensitivity studies for enhanced slant range visibility were conducted using “Empirical 
Orthogonal Functions (EOF) Regression” or “Principle Components (PC) Regression” to 
correct the effect of absorbing atmospheric constituents such as fog and aerosols. 

4.1 Clear Air Turbulence (CAT) 

Radiometric techniques for turbulence detection rely on atmospheric temperature and water 
vapor structures that depart from the background atmosphere. The radiometric techniques 
are indirect in the sense that they do not measure turbulent motion. Rather, they measure 
the temperature structure and water vapor structure of the atmosphere ahead of the aircraft 
and attempt to relate it to the likelihood that turbulence will be experienced. 

Researchers have been studying passive IR detection of turbulence for decades. For 
example, in the 1970’s, Peter Kuhn developed an IR radiometer that was flown and tested 
on three NASA aircraft in the 1979 Clear Air Turbulence Flight Test Program [7]. Kuhn’s 
passive IR system detected CAT by detecting water vapor variability ahead of the aircraft. 
(In non-CAT conditions, the water vapor content is fairly constant.) Kuhn claimed that the 
detector had a false-alarm rate of less than 8% in clear air at altitudes above 4.4 km with 
a range of 60 km. However, the system was not tested with clouds in the vicinity [7]. 

In the study reported here, FLI simulations were conducted using three-dimensional 
temperature and water vapor fields associated with CAT. Techniques for determining the 
distance and warning time to an aircraft encounter with CAT were defined. The spectral 
resolution, range, and instrument noise requirements were defined for an FLI to be used to 
detect CAT and to determine its distance from the aircraft (i.e., warning time) by both the 
temperature and water vapor radiance signal techniques. 

A high spectral resolution spherical shell radiative transfer model, based on Line by Line 
Radiative Transfer Model (LBLRTM) molecular optical depths, was developed to study 
phenomena at some distance from the aircraft. A pictorial display of the path seen by an 
FLI aboard an aircraft is shown in figure 1. 

The model takes the analytical form 

R(u,A,h 0 ) = £ OA B[u,T(a)] * [-^^-]ds (1) 

where R is radiance, TO A is top of atmosphere, v is wavenumber, A is the atmospheric 
state consisting of pressure p, temperature T, absorbing constituents (i.e., H 2 0, C0 2 , O3, 
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Figure 1. Path viewed by a forward-looking optical instrument aboard an aircraft. 

N 2 0, CO, CH 4 , 0 2 , NO, S0 2 , N0 2 , NH 3 , HN0 3 , clouds, and aerosols), h a is the height of 
the aircraft, and B is the Planck function. 

The transmittance of the atmosphere betweeen the aircraft and any position along the 
viewed path s is 



( 2 ) 


Jo 

where cr is the optical depth. The coordinate along the viewed path, s, is 


s(h) = V(r o + /0 2 -(r e + fc o ) 2 


(3) 


where h is height and r e is the radius of the Earth. 

Figure 2 shows the results of applying this model to the US Standard Atmosphere to 
simulate what an FLI would see in clear air when flying on an aircraft at three different 
altitudes. The spectral radiance is plotted versus wavenumber which would be observed by 
an FLI for aircraft altitudes of 1.5, 4.5, and 9.5 km. A Norton-Beer Weak Apodization was 
applied to interferograms simulated to a maximum optical path difference of 1.0 cm (i.e., 
an unapodized spectral resolution of 0.5 cm -1 ). 

The spectrum corresponds to a near-isothermal spectrum (i.e., a Planck function) for 
the relatively opaque 15pm C0 2 band, 650 - 750 cm -1 . A near-Planck radiance distribution 
is also observed for the strong water vapor absorption region, beyond a lower wave number 
limit ranging from 1250 - 1500 cm -1 , dependent upon instrument altitude. The strong 
emission by ozone within the 1000 - 1100 cm -1 region can also be seen in figure 2, with the 
relative amplitude of this emission band increasing with altitude as a result of the increase in 
ozone concentration and the rapid decrease of the water vapor concentration with altitude. 
The strong absorption/emission by methane and nitrous oxide near 1305 cm -1 can also 
be seen from the fact that the spectral radiance distribution approaches a Planck radiance 
spectrum as the altitude of the sensor decreases. 

Figures 3, 4, and 5 show plots of the relative contribution from the atmosphere (assuming 
US Standard Atmospheric conditions) as a function of distance from the aircraft. Each plot 
indicates the total radiance at three different altitudes, 1.5, 4.5, and 9.5 km. Contributions 
greater than 100% are a consequence of ringing in the spectrum that remains after the 
weak apodization process that was performed. The higher the aircraft, the greater the 
contribution from spherical shells at greater distances from the aircraft. For the sensing of 
turbulence from the 15 pm C0 2 band, most of the signal (> 60%) is confined to the first 
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5 km shell ahead of the aircraft for wave numbers smaller than 722 cm -1 , decreasing with 
increasing wave number. In the “window” region of the spectrum, 800 - 1200 cm -1 , the FLI 
observes contributions from the atmosphere at distances greater than 80 km when flying at 
high altitude (e.g., 9.5 km). 

This long range detection capability makes the technique useful for the detection of 
hazardous aerosol (e.g., volcanic ash clouds) contributions, providing avoidance warning 
times in excess of 5 minutes for a fast flying commercial jet aircraft. At middle and low 
atmospheric altitudes (e.g., 4.5 km and 1.5 km) the visibility of clouds and aerosol layers 
with distance decreases, being greatest for a narrow spectral region near 1000 cm -1 and the 
1100 - 1150 cm -1 region of the spectrum. In addition, high spectral resolution is needed to 
optimize the FLI instrument sensitivity to contributions of radiance far from the aircraft, 
which results from being able to sense the radiation between the water vapor absorption 
lines within the “window” region of the spectrum. 

4.2 Temperature Turbulence Measurement Concept 

Detecting temperature fluctuations associated with CAT can be accomplished using either a 
fully- or partially-scanned infrared interferometer. In the partial scan approach, by observing 
the CO2 resonance region of the interferogram produced by the near uniform spacing of CO2 
lines in the 13-15 /jm region, discontinuities of temperature will produce a signal within 
this narrow CO2 resonance region. The signal’s amplitude is dependent upon the amplitude 
of the non-uniformity of atmospheric temperature. The distance from the aircraft at which 
these discontinuities exist should be indicated by the spectral wavelength sub-region of the 
13 - 15 pm band where one begins to see the signal. For example, the stronger the signal 
gets near the more opaque center of the CO2 absorption band, the closer to the aircraft the 
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Figure 3. Contribution to the total radiance at 1.5 km. 
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Figure 4. Contribution to the total radiance at 4.5 km. 
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Figure 5. Contribution to the total radiance at 9.5 km. 


temperature fluctuation must be occurring. The radiometric signal will be small, but FTS 
systems are very sensitive to very small spectral features of the radiation being detected - 
in this case, the relative amplitude of the CO2 lines, which is zero for a completely uniform 
temperature scene. 

The spherical shell radiative transfer model was applied to a one-dimensional case of 
turbulence-related temperature variations observed from a commercial jet aircraft, which 
encountered a narrow zone of severe turbulence. The top plot of figure 6 shows the radi¬ 
ance emission lines associated with the temperature perturbations from a severe velocity 
turbulence encounter by a commercial passenger jet aircraft near Houston, TX, on August 
6, 2003 (shown in the time series bottom plot of figure 6). The relatively uniform spacing of 
the CO2 emission lines is evident in the spectrum of turbulent radiance perturbation from 
a horizontally uniform temperature condition. 

An analysis of the CO2 resonance in the interferogram produced by a turbulent temper¬ 
ature fluctuation during the same turbulence encounter by the commercial aircraft at an 
altitude of 9.5 km is shown in figure 7. In this figure, the top plot is brightness temperature 
difference versus wavenumber while the bottom plot is the interferogram ratio (/ — Ig)/Ig. 
The turbulence signal can be seen in the raw interferogram measurement from the ampli¬ 
tude of the CO2 resonance region produced by the relative uniform spacing of the CO2 
emission spectral radiance variations. Thus, the temperature fluctuation signal produced 
by wind turbulence can be detected in either the raw interferogram or the derived spectrum 
of radiance observed with an FLI instrument. 

Three case studies have been performed to determine both the interferogram and spectral 
radiance signal resulting from temperature fluctuations associated with clear air turbulence. 
One case is the actual temperature measurement fluctuation observed from a commercial 
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jet aircraft during an encounter with severe clear air turbulence, as described above, with 
the other two resulting from von Kdrmdn simulations for moderate to weak wind turbulence 
conditions. All of the simulations were conducted for the aircraft within a 40, 20, 10, 5, and 
0 km distance from the relatively narrow turbulent region. Figures 8, 9, and 10 show the 
three cases where the assumed aircraft altitude is 9.5 km. The distance to the turbulent 
region is associated with the wavenumber where significant spectral structure begins, moving 
towards larger wavenumbers with increasing distance from the turbulent regime. 

As shown in figures 11, 12, and 13, there is a relatively large spectral signal when the 
aircraft is close to the turbulence but this signal decreases rapidly as the aircraft becomes 
more distant from the turbulent region. Thus, the ability of the spectral, or partial interfero- 
gram amplitude, to provide significant warning time before severe turbulence is encountered 
depends on the interferogram signal-to-noise ratio that can be achieved with the FLI. 

Finally, figure 14 shows the altitude dependence on the spectral radiance and inter¬ 
ferogram produced from turbulent temperature fluctuations 20 km ahead of the aircraft 
simulated for von Kdrmdn case #2. In this situation, the typical commerical aircraft would 
have about 80 seconds of warning time. The signal maintains a nearly constant strength 
when the aircraft descends to 4.5 km but then decreases with a further descent to 1.5 km. 
The enhanced response of the FLI to turbulent temperature fluctuations at high altitudes 
is advantageous, since most severe clear air turbulence occurs at the higher altitudes. 

4.3 Water Vapor Turbulence Measurement Concept 

It may be easier to detect clear air turbulence at a greater distance using water vapor 
radiance measurements rather than CO 2 temperature measurements because the atmosphere 
may be relatively transparent due to low water vapor concentrations between the aircraft 
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Figure 7. C0 2 resonance in the interferogram. 
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Figure 8. Temperature fluctuation and brightness temperature signal. 
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Figure 9. von Kdrman case #1 simulation. 
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Figure 10. von Karman case #2 simulation. 
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Figure 12. von Kdrmdn case #1 severe turbulence regime. 
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Figure 13. von Kdrmdn case #2 severe turbulence regime. 
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Figure 14. Variation of spectral radiance and interferogram signal. 
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and the turbulent region. In the turbulent region, water vapor molecules are injected into 
the otherwise relatively dry upper troposphere, creating a water vapor anomaly that can be 
detected by the FLI in the water-vapor radiant emission spectral regions. 

In figure 15, the flight track of the NOAA G-IV Research Aircraft on 6 December 2003 
is shown on the left side. This picture shows the location where severe turbulence was 
encountered by the NOAA G-IV research aircraft flying in clear air over convection, over 
the Atlantic ocean off the northeast coast of the US. The image on the right side is false 
color temperature and moisture data from a satellite depicting the conditions during the 
flight. The bright red and orange indicate severe convection. 



Figure 15. NOAA G-IV aircraft flight track on 6 December 2003. 

Figure 16 is a series of plots of flight data from the NOAA G-IV flight on 6 December 
2003. Each plot is a time series of different parameters measured from the research aircraft. 
The first three plots are the u (east-west), v (north-south), and w (vertical) components 
of wind, followed by temperature, true air speed, and altitude, and finally aircraft heading 
and relative humidity. The plots indicate that there is a noticeable fluctuation in relative 
humidity at the point of the severe turbulence encounter. As a result, there should also be 
a noticeable signal in FLI radiance spectra due to this severe turbulence signature in the 
horizontal water vapor distribution. 

Figure 17 shows the spectra that would be observed for the G-IV encounter with an FLI 
observing the radiance spectrum from 680 - 1600 cm -1 . The water vapor signal of turbulence 
in the 1300 - 1600 cm -1 region is very strong (± 2 K) as early as At = 0, which is about 
150 seconds before the severe turbulence encounter. Moreover, this signal seems to persist 
over a relatively large spectral region. This persistence indicates that the spectral resolution 
of the FLI will not need to be as high for detecting the water vapor signal associated with 
severe turbulence as that required to detect the temperature signal of severe turbulence. 
The spectral resolution required for sensing the water vapor signal of turbulence is < 2.5 
cm -1 , whereas the spectral resolution required for sensing the temperature signal is < 0.625 


4.4 Volcanic Ash 

Approximately 60 volcanoes worldwide are active during a typical year, and in the North 
Pacific Region of the United States alone, volcanic ash is present an average of 4 days 
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Figure 16. Time series data from the NOAA G-IV severe turbulence encounter. 
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a year where most jet aircraft fly [8]. Volcanic ash can cause jet turbine engine failure. 
During the time period from 1973 until 2000, about 100 reported incidents occurred where 
aircraft encountered volcanic ash [2]. The Geophysical Institute at the University of Alaska 
Fairbanks recently completed an extensive study of volcanic activity in the North Pacific. 
In this study, 100 years of volcanic eruption data was used to run thousands of PUFF model 
simulations of eruptions of volcanoes in the area [9]. 1 The result was a reliable set of average 
data describing trends in ash migration from Northern Pacific volcanoes. They found that 
only 35% of volcanic eruptions in that area over the past 100 years caused ash to reach 
above 8 km or about 26,250 feet [10]. 

An Australian company, Tenix Pty Ltd, holds the patent on an airborne passive IR 
radiometer originally developed by Dr. Fred Prata when he was employed by the Com¬ 
monwealth Scientific and Industrial Research Organisation (CSIRO) located in Melbourne, 
Australia. 2 The system, named the Airborne Hazards Detection System (AHDS) uses a 
two-band detection technique similar to satellite detection techniques for volcanic ash. 

In the study reported here, investigations have been conducted into means to detect 
the distance and warning time for an aircraft encounter with a volcanic ash cloud using 
window region spectral radiance emission observations. The spectral resolution, range, and 
instrument noise requirements have been defined for an FLI to be used to detect a volcanic 
ash cloud and to determine its distance from the aircraft (i.e., warning time). Two different 
algorithms for determining the distance to an aerosol cloud were examined. The algorithms 
are the “CO2 Slicing” technique and the ratio of “window” cloud radiance relative to the 
clear air. 


4.4.1 “CO2 Slicing” Technique 

The “CO2 Slicing” technique expresses the distance to an object in terms of the optical 
depth of CO2 between the sensor and the object. CO2 is a uniformly mixed gas so the 
optical depth is proportional to the distance between the sensor and the emitting object. 
The optical depth, /, and therefore distance, is determined from equation 4. 


/ 


j^doud _ j^clear 

j^cloud _ j^clear 


(4) 


where is cloud affected CO2 channel radiance, R$ ear is the clear air CO2 channel 

radiance, Rf oud is cloud affected window channel radiance, Rf ear is the window channel 
radiance. 

In this preliminary study, the ash cloud was assumed to have a thickness of 100 km and 
an optical depth of unity. The distance predictor (ratio of CO2 channel radiance to window 
channel radiance) was produced for cloud distances of 0, 5, 10, 20, 40, and 80 km from the 
aircraft FLI. Figure 18 is a plot of the ratio of carbon dioxide to window channel radiance, 
relative to clear air, as a function of cloud distance for aircraft altitudes of 1.5, 4.5, and 9.5 
km. The plot shows the result where the 791.8 cm -1 Q-branch is used as the CO2 channel 
and 960 - 964 cm -1 is used as the “window” reference. There is little sensitivity to cloud 
distance except at 1.5 km and for the unapodized (i.e., full spectral resolution) radiance. 
This result implies that the CO2 channel used for the process is too transparent (i.e., optical 
depth of CO2 between the instrument and the cloud is too close to zero) to be useful for the 
cloud distance determination. As a consequence, more opaque CO2 regions should be used 
to determine the distance between the FLI and the ash cloud. 

The detection of aerosols (such as volcanic ash) with a moderate spectral resolution 
interferometer is made possible by a unique spectral signature that an aerosol produces 


1 PUFF is a model that simulates the movement of airborne volcanic ash in near real-time when an 
eruption has occurred. 

2 Dr. Prata is currently employed at the Norwegian Institute for Air Research. 
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Figure 18. The ratio of carbon dioxide to window channel radiance. 


in the window region relative to the spectral signature due to water droplets or ice crystal 
clouds, shown in Figure 19. This figure shows radiance spectra within the 11-12 /j ,m window 
region for a downward viewing spectrometer calculated for two different aerosol particle sizes. 
The aerosol signature is the negative slope in the radiance distribution between 11 and 11.5 
m where the absorption decreases with increasing wavelength. For water and ice clouds, 
and for water vapor, the absorption increases significantly with wavelength in this same 
spectral region. Thus, the existence of aerosols should be easily detected from a spectral 
radiance emission spectrum as would be observed with an FLI. 

The optical depth spectrum of volcanic ash is shown in Figure 20, which has a total 
optical depth of unity at 1100 cm -1 . The radiance and brightness temperature spectra that 
the FLI would observe at 1.5 km, 4.5 km, and 9.5 km with and without such an ash cloud 
are shown in Figure 21, assuming the aircraft is at a distance of 100 km from the ash cloud. 
It can be seen that the volcanic ash absorption alters the spectral distribution of radiance 
quite dramatically in the 800 - 1200 cm -1 spectral region, making the existence of a volcanic 
ash cloud easily detectable when the aircraft is flying above two kilometers in altitude. Near 
the surface (e.g., at the 1.5 km altitude), there is strong absorption by water vapor that 
tends to mask the aerosol spectral signature. In principle, the absolute magnitude of the 
radiance observed within the 800 - 1200 cm -1 spectral region should be proportional to the 
concentration of atmospheric aerosol ahead of the aircraft, whereas the amplitude of the CO2 
lines near 795 cm -1 , as well as the relative amplitude of the water vapor lines throughout 
the window region, should be indicative of how distant the volcanic ash cloud is from the 
aircraft. That is, the closer the aerosol layer is to the aircraft, the weaker the amplitude of 
the CO2 and water vapor emission spikes seen in the spectral distribution of radiance. The 
absolute distance can probably be calculated using the same CO2 slicing technique used to 
determine cloud altitude from satellite spectral radiance measurements [11]. 
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Figure 20. Volcanic ash optical depth spectrum. 


26 








































Radiance Spectra (with /without volcanic ash cloud) 



Brightness Temperatures (with / without aerosol cloud ) 



Figure 21. Brightness temperatures (with / without volcanic ash cloud). 
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Figure 22 shows the difference between the radiant brightness temperature spectra ob¬ 
served with and without a volcanic ash cloud atmospheric brightness temperature minus the 
clean atmospheric brightness temperature within the FLI instrument’s field of view. For 
this situation, where the aircraft is 100 km distant from a relatively thin (i.e., the maximum 
spectral optical depth is unity) volcanic ash cloud there is still a very large volcanic aerosol 
signal in the FLI brightness temperature difference spectrum, being in excess of 50 K near 
10.5 and 11.5 /un. Notice that the relative amplitude near the wavenumbers of 950 cm -1 
and 1100 cm -1 reverses between high and low aircraft altitudes. This reversal in brightness 
temperature difference is due to the strong decrease in water vapor emission with increas¬ 
ing altitude at 950 cm -1 , relative to 1100 cm -1 , where it is relatively small and therefore 
much less dependent upon altitude. At 1100 cm -1 , the brightness temperature dependence 
with altitude is more closely related to the volcanic ash emission profile, or volcanic aerosol 
concentration, since the volcanic ash optical depth is a maximum at this wavenumber (see 
Figure 20). There may be valuable information on the distance of the aircraft from the 
aerosol layer in the relative amplitudes in the brightness temperature signals at 950 cm -1 
and 1100 cm -1 , since the water vapor emission at 950 cm -1 should increase, relative to 
1100 cm -1 , with increasing distance from the aerosol layer. 



Wavenumber, (cm' 1 ) 


Figure 22. Difference in radiant brightness temperatures. 

Figure 23 show differences in brightness temperature spectra when the ash cloud is 
located at distance L from the aircraft (L = 5, 10, 20, 40, and 60 km) in comparison with 
the corresponding brightness spectrum when L = 0 km. Ash cloud parameters: thickness 
along the viewing path is 20 km, maximum optical depth for the whole cloud is 0.2 at A = 9 
im i. The total amplitude of the brightness temperature differences here and their large-scale 
shape are clear indicators of the ash cloud presence ahead of an aircraft. They can be used 
to retrieve the optical depth of the cloud (i.e., concentration of aerosol particles). At the 
same time, the amplitudes of the H2O, CO2, and O3 lines, visible in 750 - 1300 cm -1 range, 
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carry information about distance to the cloud, because they depend on the total air mass 
between an aircraft and an ash cloud. The farther the ash cloud is from the aircraft, the 
greater the amplitude of the molecular spectral lines. 



4.4.2 “Window” Channel Technique 

The second method for determining ash cloud distance involves simply using the ratio of 
“window” cloud radiance relative to the clear air radiance to determine the ash cloud dis¬ 
tance from the aircraft at altitudes of 1.5, 4.5, and 9.5 km. The technique is based on the 
fact that, because of the curvature of the atmosphere, the altitude of an object viewed along 
a straight-line path will vary with its distance from the instrument. Also, an ash cloud 
observed along the viewed path will raise the radiance that would otherwise be observed for 
a cloud-free path to the cold space background, since the ash cloud radiance is a function 
of both its temperature and its optical depth. Thus, assuming an opaque cloud, the radi¬ 
ance will be a function of its altitude, and therefore distance from the aircraft instrument, 
because of the atmospheric temperature dependence on altitude. 

Figure 24 is a series of plots of the “window” cloud radiance relative to the clear air 
radiance as a function of aerosol distance from the aircraft. Equation 5, is used to compute 
the cloud radiance ratio, which shows a very strong distance signal in this measurement. 
To implement this technique, however, a table of this window channel ratio as a function of 
cloud distance for different standard atmospheric conditions and different cloud altitudes is 
required. A correction for the deviation of the actual atmospheric condition from the stan¬ 
dard atmosphere could be estimated from the actual temperature and humidity measured 
at the aircraft. The process would be similar to estimating the surface pressure needed to 
obtain an altimeter setting for a particular aircraft in flight. 
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Figure 24. Ratio of cloud radiance as a function of distance. 


4.5 Low Slant Range Visibility 

One of the major hazards to safe flight is poor visibility during the landing portion of the 
aircraft flight. This condition is particularly true for commuter, business, cargo, and private 
aircraft, which use uncontrolled airports and which may not have precision approach capa¬ 
bilities such as an Instrument Landing System (ILS), Microwave Landing System (MLS), or 
Global Positioning System (GPS). At uncontrolled airports, there may be other obstacles 
(e.g., other aircraft) on the runway, which may not be seen by the pilot during a landing 
approach when Instrument Meteorological Conditions (IMC) exist due to poor visibility. 

An FLI imaging spectrometer can be used to observe the scene ahead of and below 
the aircraft from multi-spectral images of the radiance spectrum emitted by these scenes. 
Because the optical depth of the atmosphere between the scene and the aircraft produces 
unique signatures in the radiance spectrum observed for each pixel of an imaging FLI, the 
effect of the absorbing atmospheric constituents (e.g., water vapor and aerosol) on the scene 
radiance can be corrected, thereby providing the pilot with enhanced imagery of the scene 
conditions ahead of the aircraft. 

The algorithm most suitable for implementing these corrections in a real-time fashion 
is called “EOF Regression” or “PC Regression” [12], [13], [14], [15]. In this technique, the 
radiance spectrum is represented by a small number of Empirical Orthogonal Functions 
(EOFs), or Principle Components (PCs), which represent the information content of the 
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radiance spectrum in terms of a small number of statistically independent amplitudes, or 
PC scores, of the radiance spectrum representation functions (i.e., the EOFs or PCs). By 
definition, the EOFs, or PCs, are statistically optimal in the sense that no other set of 
spectral functions can be used to represent the observed radiance spectra with as small a 
number of independent variables (i.e., the EOF amplitudes or PC scores). Thus, EOF or 
PC expansions provide the most efficient representation of any particular spectral radiance 
data set. 

The algorithm to produce an enhanced, or clear, image of the background from the FLI 
imaging spectrometer data involves predicting for each detector array pixel the background 
scene temperature from the spectral radiance distribution observed by that detector element. 
The prediction equation takes the form of a linear regression equation where the predictand 
is the background surface temperature and the predictors are the amplitudes, or PC scores, 
of the first few EOFs, or PCs, used to represent the observed radiance spectrum. As will 
be seen, only three or four functions, depending on spectral resolution, are required to fit 
the window region spectrum used to infer the surface background temperature. 

The EOFs and coefficients of the prediction regression equations are obtained from 
a large statistical sample of FLI radiance spectra simulated for a variety of surface air 
temperature, surface skin temperature, aerosol, and surface visibility conditions using the 
LBLRTM described earlier in this report. As is common practice in the use of EOFs for 
remote sensing [16], the Eigenvectors are computed without the addition of noise to the 
simulated spectra. However, regression equation coefficients are derived using EOF ampli¬ 
tudes obtained from radiance spectra with random noise added at a level expected to exist 
for actual observations. A detailed mathematical description of the algorithm is provided 
below. 

The algorithm is based on Principal Component Analysis (PCA) of IR spectra obtained 
for slant viewing of the ground from an aircraft. To predict the temperature of a ground 
object, the following regression scheme is used 

T p = T 0 + f-P (6) 

where where T 0 is the surface air temperature, f is the vector of decomposing coefficients 
(PC-scores) for a given measured spectrum, R{v ), on some pre-calculated basis of Empirical 
Orthogonal Functions (EOF), {£',}, and P is a vector of predicting coefficients. 

n = R{v) ■ Ei(u) ( 7 ) 

where the subscript, i, refers to the EOF or PC number. Ei(y) and P are calculated from 
the specially-designed training set of spectra. 

LBLRTM was used to generate the basic set of spectra, R(v, To), for a given atmospheric 
model (the US standard atmosphere), the geometry of observations (altitude h and zenith 
angle a of the viewed scene), for nominal ground temperature To, and for different models 
of aerosol loading (urban and rural aerosols) with visibilities of 5, 4, 3, and 2 km (i.e., 8 
cases total). 

To recalculate the spectra for another possible ground surface temperature, T s , the 
following equation was used: 

R(v, T s ) = R{v , T 0 ) + t(v) ■ [B(u, T s ) - B(v, T 0 )] (8) 

where t(v) is total transmittance of the atmosphere for a given aerosol model, and B is the 
Planck function. 

For each of 11 basic ground temperatures selected in one degree increments between the 
nominal ground temperature, plus and minus five degrees (i.e., T s = To — 5 K, ..., To, ..., 
Tq + 5 K ), 100 Gaussian random ground temperatures (with mean = T and a = 0.5 K) 
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were generated, and then corresponding spectra were computed using equation (8). Thus, 
the full training set has 8800 spectra [he., N = 8800 = (2 aerosol models) x (4 visibilities) 
x (11 basic ) x (100 random perturbations of each ground temperature)] each spectrum 
having a different ground temperature and a varying aerosol loading. 

Figure 25 shows 16 example FLI spectra computed for a spectral resolution of 2.5 cm -1 
for a ground temperature equal to the surface air temperature of 288.2 K. Brightness temper¬ 
ature spectra, computed for the US Standard atmosphere for four different aerosol models 
(urban aerosol, rural aerosol, tropospheric extinction, and radiation fog) and four different 
visibilities (5 km, 4 km, 3 km, and 2 km) are shown in the figure. The aircraft altitude 
is 1.22 km and the instrument is viewing the ground at an angle of 20° below the horizon 
(atmospheric path between the aircraft and the surface is 3.6 km). There are four different 
extinction models (tropospheric extinction, urban aerosol, rural aerosol, and radiation fog) 
and four different surface visibilities (5, 4, 3, and 2 km) assumed for these calculations. The 
red spectrum represents the mean of the sixteen different spectra. It can be seen that the 
discrepancy between the observed brightness temperature and the ground temperature is 
related to the amplitude of the spectral structure across the window region: the larger the 
attenuation, the smaller the spectral structure when the ground temperature is equal to 
the surface air temperature. Also, different types of absorbers possess different broad-scale 
spectral shapes. Note that radiation fog, represented by the bottom four curves in the 
figure, has a distinct maximum brightness temperature (i.e., minimum extinction) near 10 
/xm (1000 cm -1 ) while the aerosol models tend to have numerous brightness temperature 
maxima across the band. 
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Figure 25. Brightness temperature spectra, at 2.5 cm 1 spectral resolution. 

Figure 26 shows the FLI brightness temperature spectra for a rural extinction aerosol 
for 11 different ground temperatures ranging between 270 K and 310 K. The surface air 
temperature is 288.2 K and the surface visibility is 3 km. The aircraft altitude is 1.22 km and 
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the view angle is 20 ° below the horizon. Figure 26 is presented to illustrate the dependence 
of the spectral variability on the discrepancy between ground temperature and the surface 
air temperature. As can be seen, the larger the discrepancy between the surface air and 
ground temperatures, the greater the spectral variability of the FLI window spectrum. Also 
note that the absorption line features flip when the sense of the difference between the 
ground and air temperature changes. It is both the large scale and small-scale spectral 
features of the FLI spectra that relate the observed brightness temperature to the actual 
ground temperature, thus enabling the correction of FLI brightness temperature imagery 
for the attenuation by the intervening atmosphere. 
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Figure 26. FLI brightness temperature spectra for a rural extinction aerosol. 

Empirical Orthogonal Functions: The training spectral radiance data set described 
above is used to calculate a mean spectrum 

RM = ^Y J R n{y, T n) (9) 

and a statistical covariance matrix 
1 N 

c°v itj = x - *>(*)! • t^i) - *0^)1 (10) 

where N is the number of radiance spectra in the training data set. A set of empirical 
orthogonal functions {Ej} is calculated as eigenvectors of the covariance matrix, cov, ordered 
in decreasing magnitude of the eigenvalues. 
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Figure 27 shows the FLI brightness temperature eigenvectors, for a spectral resolution 
of 2.5 cm -1 , computed from a sample of 4400 urban aerosol and ground temperature con¬ 
ditions, assuming the US Standard Atmosphere. The eigenvectors are computed for an 
aircraft altitude of 1.22 km and a view angle of 20 degrees below the horizon. The spectral 
resolution is 2.5 cm -1 . As can be seen the different eigenvectors portray different degrees 
of spectral variability. Thus, the amplitudes of these EOFs (i.e., the rys of equation (7), 
above), carry the information regarding the magnitude of the spectral variability of the 
observed brightness temperature spectra. 



Figure 27. FLI brightness temperature spectra for urban aerosol conditions. 


Figure 28 shows the eigenvalues, A, (the sum of which is the total variance) associated 
with the eigenvector representation of 8800 radiance spectra computed for both rural and 
urban aerosol conditions. The surface visibility ranges between 2 and 5 km and the ground 
temperature ranges between ± 5 K of the nominal surface air temperature of 288.2 K. 

The results are shown for both the 2.5 cm -1 and 0.6 cm -1 spectral resolutions considered 
in this study. Note that the square root of A is plotted on a logarithmic scale in order to 
differentiate the two curves for the higher order EOFs. There is additional information 
content for the higher spectral resolution, which increases with increasing EOF number. It 
is noted that even though nearly twice as many EOFs are needed to explain all the variance 
of the high resolution spectra, as compared to the lower resolution spectra, the absolute 
magnitude of the residual variance is relatively small for EOF numbers beyond three. 

Vector of prediction regression coefficients P: Having the basis of empirical orthogonal 
functions, the matrix C = |cj ;J | can be calculated, where c,,j are PC-scores for all spectra 
from the training data set; that is, 
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Figure 28. Eigenvalues of a matrix of the statistical covariance. 


Ci,j = Ri ■ Ej (11) 

where i refers to the training data set sample number and j refers to EOF number. 

The elements of prediction vector P = {pj } are computed by solving the system of linear 
equations Cij • Pj = T t — Tq, or in matrix form 

CP = AT. (12) 

Because equation (12) is an over-determined system of equations (i.e., more equations 
than unknowns), a least squares solution minimizes the sum of the squares of the differences 
between the actual ground temperatures and predicted ones. Thus, solving for that P for 
which (AT — CP) T (AT — CP) = min., the least squares solution is 

P = (C t C)- 1 C t AT. (13) 

The solution is conditioned with respect to instrument noise amplification by adding 
Gaussian random noise, with standard deviation, a, equal to the expected measurement 
noise, to each radiance measurement in the training data set before calculating the matrix 
C using equation (11). 

Figure 29 shows the root-mean-square error of predicted surface temperature for 8800 
spectra simulated for urban and rural aerosol conditions with surface visibilities ranging 
between 2 and 5 km and surface temperatures ranging between 283 and 293 K, for Stan¬ 
dard Atmospheric conditions with a surface air temperature of 288.2 K. The results for 
four different FLI instrumental noise levels are shown. The number of EOFs required to 
predict the background temperature is four, or less, depending on spectral resolution and 
instrumental noise. This result is important in that it should be easy to perform a simple 
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two, three, or even four eigenvector regression estimate of surface temperature for every 
pixel of an imaging interferometer in real time aboard the aircraft. It can also be seen that 
for realistic noise levels of 0.3 K, the background temperature can be estimated with an 
accuracy close to 0.5 K, which is only slightly poorer than the brightness temperature noise 
level. Scene contrasts on the order of a one half of a degree, or more, should be easily seen 
in the atmospheric-extinction-corrected FLI imagery. 



Figure 29. Root mean square error of predicted surface temperature. 

Figure 30 shows, for a spectral resolution of 2.5 cm -1 , the dependence of predicted 
ground temperature accuracy on the surface visibility condition. The dependence of the 
accuracy on surface visibility decreases with increasing number of eigenvector amplitude 
predictors used for the extinction correction process; however, the residual dependence is 
strongly a function of instrument brightness temperature measurement noise level. If the 
instrument noise level can be maintained at a level near 0.1 K, then accuracies better than 
0.3 K can be achieved, even under Low Instrument Flight Rules (LIFR) flight conditions 
(e.g., visibility < 2 km). 

FLI imagery simulations Regression prediction equations were specified to correct 
simulated FLI imagery for atmospheric extinction. Table 1 is a summary of the various ob¬ 
serving conditions used for the simulations conducted to determine the radiance eigenvectors 
and regression equation prediction coefficients. Ranges of values for some parameters are 
listed. 
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Figure 30. Dependence of predicted ground temperature accuracy. 


Table 1. FLI INSTRUMENT ATMOSPHERIC/SURFACE SIMULATION VALUES 


Parameter 

Values 

Wavenumber range 

800 - 1200 cm" 1 

Wavenumber resolution 

2.5 / 0.6 cm -1 (low/high) 

Atmospheric profiles 

US Standard Atmosphere 

Aerosol Models 

LBLRTM Urban Aerosol and Rural Aerosol 

Surface visibility 

5.0, 4.0, 3.0, and 2 km 

Surface air temperature (To) 

288.2 K 

Ground temperature (T p ) 

283.2 to 293.2 K in increments of 

1 ±e with a(e) = 0.5 K 

Aircraft Altitudes (h) 

h = 4000 ft. (~1.22 km) and h = 1000 ft. (~0.3 km) 

Zenith angle of observations 

a = 110° (—20° to horizon) 

Slant range to the ground target 

~ 3.6 km for h = 4000ft.; ~ 0.9 km for h = 1000ft. 

Number of predictors 

4 

Brightness temperature noise 

cr(e) = 0.1, 0.3, and 1.0 K 
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Figures 31 through 35 show imagery created from the FLI’s most transparent channel 
brightness temperature and EOF regression predicted background temperature for a variety 
of atmospheric visibility, instrument noise level, and aircraft altitude conditions. The scene 
is composed of three radiating target temperatures, grass, a concrete runway, and a vehicle 
on the runway, their temperatures indicated in the figure legends. These radiating tempera¬ 
tures might represent a night-time condition where the runway is radiatively cooler than the 
surrounding grass. The radiating temperature difference might result from the difference in 
the infrared emissivity of the concrete (or asphalt) and/or the more rapid infrared cooling of 
the runway surface at night than the surrounding grass-covered dirt. The opposite temper¬ 
ature differential between the runway and the surrounding turf might be expected to occur 
during the day when the surface is illuminated by the sun. As simulated, one would expect 
contrast between a vehicle, or other object, on the runway due to its different thermal and 
emissivity properties from those for the runway on which it is positioned. 

The EOF-regression-predicted background imagery does a remarkable job of enhancing 
the contrast between the runway and a possible obstruction on it. As expected, slightly 
better results are obtained at high spectral resolution, particularly for high noise measure¬ 
ment conditions. However, considering the fact that the measurement noise level of a 2.5 
cm -1 FLI should be smaller than that for a 0.6 cm -1 resolution FLI (a factor of ~2 for the 
same Michelson mirror scan rate) and the fact that the refresh rate would be improved by 
a factor of 4 as a result of the shorter interferogram length, the 2.5 cm -1 FLI system is pre¬ 
ferred for this particular application. Further investigations into the dependence of the FLI 
slant range imagery on spectral resolution might reveal that even poorer spectral resolution 
(shorter interferogram scan time) measurements would be sufficient, thereby enabling even 
more continuous imagery than is possible with the higher spectral resolution system. The 
EOF regression algorithm works remarkably well and has the ability to retrieve very clear 
scene condition imagery under very low visibility conditions, even over large slant ranges. 
This result indicates that an FLI imaging spectrometer would be extremely useful in the 
cockpit to assist the pilot in seeing hazardous obstructions along the flight path and for 
enabling visual landings under low visibility conditions. 

Figure 31 is a comparison between most transparent channel brightness temperature 
imagery with EOF regression enhanced imagery for an aircraft altitude of 4000 ft (slant 
range equal to 3.6 km), and instrument noise level of 0.1 K, and surface visibility of 1 km 
(~ 0.6 mi) for 0.6 cm -1 and 2.5 cm -1 spectral resolutions. The scene is composed of three 
radiating target temperatures, grass at a temperature of 288 K, a concrete runway at a 
temperature of 286 K, and a vehicle on the runway with a temperature of 290 K. 

Figure 32 is a comparison between most transparent channel brightness temperature 
imagery with EOF regression enhanced imagery for an aircraft altitude of 1000 ft (slant 
range equal to 0.8 km), and instrument noise level of 0.3 K, and surface visibility of 0.25 
km (~ 0.16 mi) for 0.6 cm -1 and 2.5 cm -1 spectral resolutions. The scene is composed of 
three radiating target temperatures, grass at a temperature of 288 K, a concrete runway at 
a temperature of 287 K, and a vehicle on the runway with a temperature of 289 K. 

Figure 33 is a comparison between most transparent channel brightness temperature 
imagery with EOF regression enhanced imagery for an aircraft altitude of 1000 ft (slant 
range equal to 0.8 km), and instrument noise level of 1.0 K, and surface visibility of 0.25 
km (~ 0.16 mi) for 0.6 cm -1 and 2.5 cm -1 spectral resolutions. The scene is composed of 
three radiating target temperatures, grass at a temperature of 288 K, a concrete runway at 
a temperature of 286 K, and a vehicle on the runway with a temperature of 290 K. 

Shown in figure 34 is a comparison between most transparent channel brightness tem¬ 
perature imagery with EOF regression enhanced imagery for an aircraft altitude of 1000 ft 
(slant range equal to 0.8 km), and instrument noise level of 1.0 K, and surface visibility of 
0.5 km (~ 0.32 mi) for 0.6 cm -1 and 2.5 cm -1 spectral resolutions. The scene is composed 
of three radiating target temperatures, grass at a temperature of 288 K, a concrete runway 
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Figure 31. Visibility = 1 km, slant range = 3.6 km, noise level = 0.1 K. 
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Figure 32. Visibility = 0.25 km, slant range = 0.8 km, noise level = 0.3 K. 
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Figure 33. Visibility = 0.25 km, slant range = 0.8 km, noise level = 1.0 K. 
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at a temperature of 287 K, and a vehicle on the runway with a temperature of 289 K. 






Figure 35. FLI settings from figure 33 but targets at higher temperatures. 


Figure 35 is a comparison between most transparent channel brightness temperature 
imagery with EOF regression enhanced imagery for an aircraft altitude of 1000 ft (slant 
range equal to 0.8 km), and instrument noise level of 1.0 K, and surface visibility of 0.25 
km (~ 0.16 mi) for 0.6 cm -1 and 2.5 cm -1 spectral resolutions. The scene is composed of 
three radiating target temperatures, grass at a temperature of 288 K, a concrete runway at 
a temperature of 287 K, and a vehicle on the runway with a temperature of 289 K. 

Based on these sensitivity studies, an imaging FLI spectrometer is very useful for pro¬ 
viding enhanced vision in the cockpit of an aircraft. The imagery can be greatly enhanced 
using a multi-spectral processing algorithm which retrieves the background scene temper¬ 
ature distribution from the observed radiance spectra. The algorithm is a simple linear 
operator and therefore could be implemented for the provision of real-time cockpit displays 
of enhanced vision of the scene ahead of the aircraft. 

The algorithm was demonstrated under a restricted set of aircraft altitude and atmo¬ 
spheric conditions. Much more work needs to be performed to demonstrate the robustness of 
the algorithm and to define the actual measurement and display characteristics of a practical 
image enhancement system. 


43 














5 Characterization Studies of Other Hazards 


Characterization studies have been performed to determine the possibility of detection of 
wake vortices, dry wind shear, and icing by an FLI instrument. Hazard phenomenology has 
been investigated that supports detection by an FLI sensor, such as temperature and water 
vapor content that were crucial to the sensitivity studies. In addition, other phenomenology, 
such as exhaust gas constituents and the optical properties of liquid water and ice, were 
examined as a further means of detection with an FLI. 

5.1 Wake Vortex 

The possibility of detecting wake vortices with an imaging FLI instrument has been inves¬ 
tigated. A study of the temperature structure, time evolution, and exhaust entrainment in 
wake vortices was conducted so that the likelihood of detection with an FLI system could 
be determined. To prevent hazardous wake turbulence encounters between aircraft, airports 
operate with the Instrument Flight Rules (IFR) longitudinal spacing requirement. The spac¬ 
ing requirement is based on factors such as the geometry of the airspace and runways and the 
size of the aircraft involved. Through the Aircraft Wake VOrtex Spacing System (AVOSS) 
project, NASA set the goal of demonstrating a technology that could provide aircraft with 
real-time wake avoidance information, thereby allowing a greater capacity of aircraft land¬ 
ings by reducing the longitudinal spacing requirement. AVOSS focused on vortices created 
during takeoff and landing [17]. 

There are several models available that describe the vorticity, shape, size and tempera¬ 
ture profiles of wake vortices created by aircraft. Temperature changes and water concentra¬ 
tions are based on the properties of fluid dynamics. Because of the interest in condensation 
trail (contrail) cloud formation and the effects on global climate, most of the data on exhaust 
entrainment is taken at or modeled for cruise altitudes. However, much of the information 
that is gathered during the study of contrails is relevant to wake vortices near the ground 
as well. 

There are four phases of contrail evolution which include Jet, Vortex, Dissipation, and 
Diffusion. The Jet stage lasts 1-20 s; the Vortex phase lasts 20-100 s; the Dissipation phase 
lasts 100-1000 s; and the Diffusion phase can last hours. 3 During the Jet phase, aircraft 
engine exhaust material is entrained in the vortices. During the Vortex phase, the vortices 
descend due to mutual induction [20], [18]. This study focuses on the Jet and Vortex phases 
as the two most relevant for vortices that affect runway spacing requirements, though some 
vortices can last several minutes in the runway environment. 

When a vortex is formed, the air in the core expands so that there is a pressure differ¬ 
ential across the swirling air. The pressure differential provides the centripetal force that is 
required to keep the air in circular motion. Therefore, the air in the core is cooled by an ini¬ 
tial adiabatic expansion. The temperature at the core of a wake vortex is typically 6 — 7°C 
cooler than the ambient air. However, this temperature differential depends strongly on the 
core size [21]. 

There have been several studies conducted that measured temperature drops in the 
vortices of an aircraft. Using anemometers, Hallock et al. often measured a 1 °C decrease in 
temperature in vortices in a study at Frankfurt/Main International Airport [22]. Adamson 
and Morrison used a scanning multi-band radiometer in 1991 to make low-resolution images 
of wake vortices, using a two-band differential signal processing technique to maximize the 

3 There is some disagreement here on whether there axe three or four phases, and what the timescales 
are. We have chosen to use the definitions outlined by Paugam et al. (2006) [18] but the definitions in 
reference [19] are similar. They define three regions of interest: jet, vortex, and plume dispersion. The time 
scales associated with these regimes are 0-10 s for the jet, 10-100 s for the vortex, and 100 s to tens of hours 
for plume dispersion. 
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system’s response at a desired range. The vortex signature was normally a cold infrared wake 
at —2 to —3°C below ambient, although the DC-10 produced a temperature difference of 
—7°C on one occasion [23]. It should be noted that both Hallock and Adamson measured 
one incident where the vortex registered a warmer core, rather than a cooler core. As 
another example, Baumgardner et al. took measurements by penetrating the contrails of 
a Lear 35 during seven flight tests to evaluate the temperature, water vapor content, and 
velocity structure of the contrails and compare the measurements to the Aerodyne wake 
vortex model. 4 They found that the increase in temperature at 500 m from the aircraft 
exhaust release point is 0.6 °C, and the increase in water vapor content was 8 ppmv. The 
measurements verified their temperature models, but showed that their models produced 
higher than actual values for the water vapor content [24]. Even with these studies, a lack 
of data still exists due to the variance in the size of aircraft and meteorological conditions. 

As the vortex ages, it weakens and decays. Vortex decay has been modeled by several 
researchers, including Greene, Sarpkaya, and Holzapfel, with the latest model accounting 
for conditions such as wind, turbulence, and proximity to the ground [25], [26], [27]. During 
decay, the core size stays the same, which means that the core pressure increases, caus¬ 
ing adiabatic heating, so the temperature differential becomes smaller. In addition, the 
whole two-vortex system sinks toward the ground as time goes on. This sinking also causes 
adiabatic heating [21]. At ground level, vortex motion and decay are very dependent on 
meteorological conditions, as well as the aircraft type and position. 

Figure 36 is an image of a vortex created by an aircraft that can be seen due to a smoke 
injection system [28]. It is apparent from the image that there is a central core of the 
vortex that is surrounded by swirling air. Aircraft engine exhaust gasses are entrained in 
the circulating air, but not in the core of the vortex. Contrails also show this phenomenon. 
The concentrations of exhaust particles in the vortex will be higher than that of the ambient 
air [19]. 

The fraction of the exhaust that gets entrained depends on whether the engines are on 
the wings or on the aircraft body. F. Garnier et al. analyzed the entrainment processes 
of aircraft engine exhaust within the wake vortex generated by the aircraft using modeling 
and simulation of a two-engine Airbus 330 and a four-engine Boeing 747. The investigation 
focused on a thirty-second timescale including the Jet and part of the Vortex phases of the 
wake vortex. The investigation found that entrainment rate is affected by many parameters 
including the jet engine location, which affects the early wake in the Jet stage, and the 
buoyancy effect in the Vortex phase after entrainment has occurred. Once entrained, the 
exhaust concentration is variable within the contrail due to stretching and distortion, as 
well as the descent of the vortex [20]. Even at the Dissipation phase, particle concentrations 
are still highly elevated over background levels of ambient species, but they are now at 
temperature and pressure conditions of the background atmosphere [19]. 

Aircraft engine exhaust consists of several species that have identifiable signatures. A 
report published by the US Environmental Protection Agency (EPA) in 1999, which eval¬ 
uated the impact of ground-level aircraft emissions generated during takeoff and landing 
on ten US cities, found that volatile organic compounds (VOCs), CO, NO x , particulate 
matter (PM), and SO2 are the five most significant species emitted from commercial jet 
aircraft that affect local air quality. During idling and taxiing when engines are operating 
at low power, VOCs and CO emission rates are highest. The highest rates of NO x and PM 
emissions occur during takeoff and climb-out. NO x emissions rise with an increase in the 
operational power of the engine, and PM emissions rise due to the incomplete combustion of 
fuel during those same operational periods. During general combustion of fuel that contains 
sulfur, S0 2 emissions are created [29]. Another publication lists major species as H 2 0, C0 2 , 
CO, and NO x and minor species as H 2 S04, and HNO3. [20]. 

4 This study was funded by the NASA High Speed Research Program. 
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Figure 36. Example of red-colored smoke entrained in wake of aircraft. 


Note that the Civil Aviation Authority in the United Kingdom hosts a website that 
provides access to the Council of the International Civil Aviation Organization’s (ICAO) 
Engine Exhaust Emissions Databank of aircraft engines that have entered production. En¬ 
gines are sorted by engine manufacturers such as Pratt & Whitney Aircraft Group and 
General Electric, and each manufacturer has a listing and corresponding datasheet of the 
various engines in production. The emissions quantified are unburned hydrocarbons (HC), 
CO, and NO x [30]. 

5.2 Dry Wind Shear 

In 1986, NASA and the FAA jointly set a timetable for developing and demonstrating a 
solution to the aviation hazard that was responsible for more than half of the US commercial 
aviation fatalities of the previous decade: wind shear. Several detection technologies were 
investigated, including radar, Doppler lidar, and infrared radiometry. During the course of 
this effort, five NASA - FAA wind shear conferences were held, with the most recent one in 
September, 1993. 

During the wind shear research effort, passive infrared (IR) detection was studied exten¬ 
sively by H. P. Adamson, who founded the Turbulence Prediction System (TPS) company to 
commercialize his infrared Advance Warning Airborne Sensors (AWAS) [28] [29]. By 1988, 
Adamson had already planned a development program including extensive in-service evalua¬ 
tion of a radiometric hazard detection system on commercial airline flights, integration with 
other avionic systems, pilot training, and FAA certification [30]. Much of Adamson’s work 
is described in the proceedings of the NASA - FAA wind shear conferences, and additional 
information is given in Appendix B of this report. 

The term dry wind shear is used here to mean low level wind shear caused by downdrafts, 
as illustrated in Figure 37 [31]. The physical basis for radiometric detection of such wind 
shear is the atmosphere’s temperature lapse rate, which causes air at higher altitude to be 
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colder. When the cold air is brought to lower altitudes by a downdraft, a large temperature 
anomaly is created. 
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Figure 37. Dry wind shear detection illustration from Adamson. 


5.3 Icing 

Ice and liquid water have very different optical properties in the infrared molecular “window” 
region of the spectrum between 10 and 12 /an. As a result, it will be possible to use the 
FLI instrument to diagnose the existence of supercooled liquid water in clouds (i.e., detect 
water droplets at sub-freezing temperatures that causes airframe icing) and to discriminate 
between ice and water on surfaces, such as runways. The temperature of the cloud can be 
observed as the brightness temperature measured with the FLI near 12 pm, where the cloud 
emissivity is a maximum, or assumed to be close to the outside air temperature observed 
independently by the aircraft measurement system. Figures 38 and 39 show that the real 
(reflection) and imaginary (absorption) parts of the index of refraction for water and ice 
differ greatly in the range of infrared wavelengths where atmospheric molecular absorption 
is a minimum (i.e., the molecular absorption “window” between 10 and 12 pm) [32], [33]. 

In particular, when viewing cloud ice crystals or an icy runway surface, with an FLI, the 
radiance observed will increase rapidly between 10.5 and 12 pm. Whereas when viewing 
water droplets or a wet surface, the radiance observed will be relatively constant between 
10.5 and 11 pm with the strong increase in radiance occurring between 11 and 12 pm. Thus, 
it is the small spectral shift in the wavelength at which the radiance emission abruptly 
increases across the atmospheric “window” region that will allow the FLI to determine 
whether cloud icing conditions, or a hazardous icy runway, may be encountered by the 
aircraft. There is a precedent for this type of discrimination; two-band IR radiometers are 
commercially available for detecting “black ice” on highways (black ice is the term for a thin 
layer of ice that is difficult to detect visually). 
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Figure 38. Real index of refraction for water and ice in infrared. 
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Figure 39. Imaginary index of refraction for water and ice in the infrared. 
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6 Instrument Design Concept and Specification 

Based on the results of the sensitivity studies and characterizations shown in the previous 
sections of this report, a table of FLI radiance measurement requirements was developed. 

Table 2 was developed to guide the development or selection of an instrument to observe 
various aviation weather hazards in future experimental studies. As discussed previously, 
there are two techniques for using FLI radiance data to detect turbulence: (1) the use 
of the C0 2 emission signal produced by the spatial variance of atmospheric temperature 
associated with the adiabatic heating and cooling resulting from turbulent updraft and 
downdrafts of the atmosphere, and (2) the use of the water vapor radiance signal produced 
by the temperature and associated humidity variations associated with turbulent updraft 
and downdrafts of the atmosphere. 

In Table 2, the instrument spectral measurement and noise requirements for the use of 
these two techniques are quite different. For the temperature turbulence using the partial 
interferogram approach, the spectral region must contain the C0 2 emission spectral region 
(690 - 715 cm -1 ) with a spectral resolution of 0.625 cm -1 , or better, in order to resolve the 
C0 2 emission line spacing of 1.5 cm -1 . As discussed earlier, Nyquist sampling of the C0 2 
emission line spacing will produce a resonance in the 0.55 to 0.75 cm optical delay region 
of the measured interferogram. The partial interferogram signal amplitude is related to the 
degree of temperature variability (i.e., related to turbulence) ahead of the aircraft. For this 
type of detection, a single field of view/single pixel spatially scanning instrument would 
suffice, but better detection would certainly be provided by an imaging Fourier Transform 
Spectrometer (FTS) with a much larger field of view and smaller pixel size than that specified 
in Table 2. 

The water vapor radiance turbulence measurement can be achieved with a lower spectral 
resolution (~2.5 cm -1 ) FTS covering the water vapor emission bands of 1200 - 1600 cm -1 
and/or 1600 - 2000 cm -1 . The reason for the potential of using one or the other of these 
water vapor bands is due to the fact that their emission characteristics are similar since they 
represent the longwave and shortwave sides of the 6.3 m band centered at a wavenumber 
of 1600 cm -1 . However, there are some subtle advantages to using one side of the band, 
versus the other, for the water vapor turbulence signal detection. The shortwave side has 
the advantage being free from absorption by spectrally overlapping trace gases, such as 
N 2 0 and CH 4 , which overlap the H 2 0 lines along the longwave side of the band. On the 
contrary, because of the rapid fall-off of Planck radiance, for a given scene temperature, 
with increasing wavenumber (decreasing wavelength), the signal-to-noise ratio of an FTS 
will naturally be higher for measurements along the longwave side of the band, the degree 
of which is dependent upon the high wavelength cut-off of the detector used. In any case, 
the choice of the band to be used will depend upon the ultimate measurement objectives 
for the FLI, which will determine the total measurement spectral range, and the properties 
of the particular optical elements and detectors available for the development of an FLI 
instrument. 

The volcanic ash detection requirements are similar to those for water vapor radiance 
turbulence, but the noise requirement is less stringent as a result of the aerosol being a much 
stronger absorber/emitter than gaseous water vapor. The band to be used for volcanic ash 
detection is extended to 1400 cm -1 in order to also observe the volcanic emissions of S0 2 , 
which is a strong absorber and emitter between 1325 and 1375 cm -1 . 

The three aviation hazards which require an imaging FTS for detection are slant range 
visibility, icing, and wake turbulence, though the detection of all other aviation hazards 
would benefit from the use of a spatially scanning imaging spectrometer, rather than a 
single pixel instrument. For slant range visibility, it has been previously shown in Section 

4.5 that effective enhancement can be achieved with a relatively low spectral resolution (i.e., 

2.5 cm -1 ) instrument operating in the thermal “window” (i.e., 800 - 1200 cm -1 ) region of 
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Table 2. FLI RADIANCE MEASUREMENT REQUIREMENTS 
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the spectrum, although improvements with higher spectral resolution (e.g., 0.6 cm -1 ) were 
shown for very low visibility conditions. 

For wake turbulence detection, the spectral range of the imaging FTS must be increased 
to obtain contributions from carbon dioxide and water vapor since it is the spatial structure 
of these emissions which delineate the temperature and water vapor spiral patterns associ¬ 
ated with the wake vortices produced by an aircraft. The field of view needs to be relatively 
large to capture any wake, which might affect an aircraft on takeoff, landing, or in-flight. 
The large field of view, coupled with the very small pixel size (i.e., high spatial resolution) 
drives the requirement towards a fairly large focal plane array (e.g., ~ 500 x 500 pixels) 
imaging device. A relatively high image refresh rate is required to enable expanded spatial 
coverage by stepping the instrument’s telescope or an associated scene scan mirror, while at 
the same time providing near continuous imagery at high spatial resolution of the obstacle 
as the aircraft moves closer to the scene (i.e., the runway during landing or a wake vortex 
during takeoff or in-flight). 

Although this measurement capability is challenging given the current state of the art, 
there are instruments which have the spectral imaging capability sufficient for experimentally 
demonstrating the concept. Also, the array size and associated field-of-view requirement 
could be relaxed by increasing the spatial stepping rate of the telescope or scene scan mirror 
of the imaging FTS. Detailed trade studies need to be performed by potential instrument 
manufacturers in order to optimize the FLI design. Such engineering studies are beyond the 
scope of the current theoretical study effort. 

7 Future Research Needs 

The sensitivity and characterization studies reported here addressed the phenomenology that 
supports detection by the FLI, in other words, which airborne hazards may be detectable, 
and how. The simulations revealed detectable signatures in both the interferogram and 
spectral domains for temperature anomalies, and detectable spectral signatures for humidity 
anomalies caused by CAT and wind shear. A literature search on wake vortex properties 
suggested that an FLI will be able to provide images of vortices. A simulation of multi¬ 
band imaging in low slant-range visibility conditions illustrated the potential of the FLI to 
provide essentially unobscured IR images of runway scenes as viewed from the air during an 
approach. In addition, the ability of an FLI to distinguish volcanic ash clouds from water 
and ice clouds was verified, and its ability to distinguish ice from liquid water on surfaces 
such as runways was investigated. Techniques for determining the range (and hence warning 
time) to hazards were also demonstrated. 

The results of these studies are promising, but they are only a first step in assessing 
the full potential of an FLI to serve as a multi-hazard airborne sensor. Investigators must 
also develop a means to estimate the intensities of the hazards posed to an aircraft; de¬ 
velop instrument requirements for a commercial hazard warning sensor 5 ; develop robust 
algorithms to relate hazard levels to sensor readings; and perform validation tests. In order 
to accomplish these goals, further research is required in four general categories including 
ground-based measurements, airborne measurements, further sensitivity studies, and simu¬ 
lations. 

7.1 Ground Measurements 

Ground-based measurements can be performed at lower cost than airborne measurements 
and may be achievable for several hazards including mountain waves (turbulence induced by 

5 Discussions have been held with two companies that have either a commercial or prototype system that 
could be used in future field tests. Information regarding these systems is included in Appendix C. 
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orography), wake vortices, low slant-range visibility, and volcanic ash. Mountain waves offer 
an excellent opportunity for ground-based measurements of turbulence, because they are 
predictable and fairly well understood. Imagery of mountain waves with an FLI may provide 
a good first comparison with the sensitivity studies. Imagery in low visibility conditions can 
be collected on any scenario where terrain or an elevated structure provides a slant viewing 
path to the earth, in a location where fog and haze occur. An airport would be the ideal 
location, and a control tower may provide the viewing platform. The airport weather data 
would document the weather conditions, and visible-light imagery would provide additional 
characterization. 

Data has been collected on volcanic ash and SO2 in the past, but no one has ever acquired 
multi-band IR imagery of a volcanic ash cloud viewed horizontally, as it would be viewed 
in flight. Such a view may be achievable with the sensor located on an adjacent mountain. 
This scenario is unpredictable because it is dependent on having an active volcano with 
an accessible adjacent mountain, but researchers should be prepared for it in case such an 
opportunity arises. 

7.2 Flight Measurements 

The number of atmospheric data sets available for the sensitivity studies reported here was 
very limited, and details of the atmospheric structures and fields associated with hazards 
remains largely unknown. For this reason, there is no substitute for future flight tests. They 
will be essential for establishing sensor performance. However, such flights must include 
truth data that can be used in the simulation environment described below to provide a 
physical understanding of the infrared signatures and to identify any issues. In order to 
keep the cost of airborne tests affordable, it may be necessary to align with larger projects. 
There are several ongoing research projects that may offer candidate airborne platforms 
along with relevant truth data, such as: 

• The NOAA G-IV is operated in a Winter Storm series during which truth data is 
acquired with drop sondes. 

• The NOAA HIAPER is sometimes flown with a wind-scanning Lidar. 

• The DLR Falcon is flown with a wind-scanning Lidar. 

• A UK Met aircraft located in Iceland is flown in turbulence studies. 

7.3 Further Sensitivity Studies 

The pilot studies reported here have successfully demonstrated that the spectral radiance 
measurements by an FLI mounted on an aircraft would be able to sense radiometric signa¬ 
tures associated with aviation weather hazards to safe flight. It has also been shown that, 
for the case of enhancing slant range visibility, Eigenvector (EOF) regression is a powerful 
method for transforming FLI spectra into useful physical products that can be used by the 
flight crew to insure safe flight. It is likewise expected that EOF regression can also be used 
to enable rapid and accurate predictions of other aviation safety hazards such as turbulence, 
wind shear, volcanic ash, wake turbulence, and icing. 

The advantage of using EOF regression is that it enables the entire spectrum of infor¬ 
mation content to be used for the sensing of these phenomena, thereby multiplexing the 
emission signals from all the atmospheric gases of interest as well as the emission signals 
from clouds and aerosols. For example, it was shown that for turbulence the signal comes 
from the temperature and water vapor contribution to the emitted radiances in two different 
portions of the emission spectrum (i.e., the CO2 emission region near 15 m and the H 2 0 
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emission region near 6.3 /an). Furthermore, for the case of turbulence, the signal is rela¬ 
tively small so that a very high signal-to-noise ratio of the measurement variable is required. 
The advantage of the use of EOF regression is that by extracting the information content 
from the entire spectrum, the signatures of turbulence inherent in the large number (i.e., 
1000 or more) of FLI spectral radiance measurements, are combined into a small number 
(i.e., 2 - 4) of EOF predictors. As a result, the signal-to-noise ratio of the FLI detection 
system is enhanced by an order of magnitude, or more, by using EOF regression to retrieve 
the geophysical characteristics of the otherwise small radiance turbulence signal. A simi¬ 
lar signal-to-noise advantage will be gained for the detection of all other weather related 
hazards as well. Thus, the important next modeling step in the FLI applications research 
is to develop and apply the EOF regression technique to predicting the intensity, distance, 
and time to encounter for turbulence, volcanic aerosol, wake turbulence, and icing. The 
technique should be exercised for a wide variety of atmospheric situations to demonstrate 
the robustness of the algorithm. 

7.4 Simulation Environment 

It is not feasible to fly a sensor in all possible conditions, so simulations generally play a 
crucial role in characterizing phenomena and developing hazard detection algorithms and 
estimates of metrics such as probability of detection and false alarm rate. In addition, all 
flight and ground test data (both sensor data and truth data) should be processed using 
the same tools as those used for simulated data. Consistency in data processing will aid 
in understanding the phenomena and the mechanisms that create aircraft hazards. It will 
also help in understanding anomalies or inconsistencies between test data and simulations. 
Knowledge gained in the comparison of test data and simulations will both help validate 
the simulations and guide the planning of future flight tests. 


8 Conclusions 

In the study reported here, remote sensing techniques developed for space-based sensors were 
applied to aviation hazard detection problems in order to assess the feasibility of using such 
techniques in airborne sensors. Sensitivity studies were performed on clear air turbulence, 
volcanic ash clouds, and slant range visibility. Characterization studies were conducted 
on wake vortices, dry wind shear, and icing to determine if future sensitivity studies were 
warranted. Sensitivity studies were conducted using a limited set of temperature, wind, 
and humidity fields. The studies were performed spectrally with high enough resolution 
to resolve individual spectral lines of molecules such as CO2, so as to support the use of 
interferometric sensors such as the FLI. The simulations were performed for single lines of 
sight from the sensor to/through the hazards. 

CAT and Dry Wind Shear. Sensitivity studies for CAT detection included two tech¬ 
niques: sensing temperature fluctuations using the 650 - 780 cm -1 C0 2 band and sensing 
water vapor fluctuations in the 1300 - 1600 cm -1 region. The simulations revealed detectable 
signatures in both the interferogram and spectral domains for the CO2 temperature tech¬ 
nique at flight altitudes, although the maximum detection range appears to be about 5 km 
(assuming an instrument resolution of 0.3 K), which is not as large as early investigators 
have claimed, but humidity anomalies in the one case for which data exists were detectable 
from 30 km, or with a warning time of approximately 2 min, assuming a lower instrument 
sensitivity of 1.0 K. Wind shear near the ground is associated with the downward motion of 
cooler air from higher altitudes, and the anomalously low air temperatures should be easily 
detectable with an FLI during a landing approach. 

Volcanic Ash. The ability of an IR imager to distinguish volcanic ash clouds from 
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water and ice clouds was verified in the 800 - 1100 cm -1 region, using an ash optical depth 
spectrum obtained from the literature. The optical properties of volcanic ash produce a 
large signature in the 8-12 /xm window region, thereby making ash detection possible with 
the FLI. 

Low Visibility. A simulation of multi-band imaging in low slant-range visibility con¬ 
ditions illustrated the potential of the FLI to provide essentially unobscured IR images of 
runway scenes as viewed from the air during an approach. This capability will give pilots 
a clear view of runways even when an obscurant such as a ground fog is present, and it 
will provide high quality images to support automated detection of airborne and ground 
obstacles or other hazards. 

Wake Vortices. Typical wake vortex properties obtained from a literature search 
included the core temperature as well as the concentrations of exhaust gas constituents 
(including water vapor) that become entrained in the vortex. Vortex core temperatures 
have been measured to be several degrees below ambient temperature and this temperature 
difference should be detectable, which suggests that an FLI might provide real-time imagery 
of the motion and evolution of wake vortices. The entrained gasses, including water vapor, 
could provide additional means of detection with an FLI. Detection of wake vortices by FLI 
techniques is possibly the most challenging application considered due to the relatively small 
size of the phenomena and the geometry of the problem which would have the FLI looking 
down at a scene with ground, runway, and other objects as background during landing, the 
phase of flight that would be of most interest. 

Icing. The optical properties of liquid water and ice are quite different in the spectral 
regions where the FLI will operate. This fact will give the FLI a capability to distinguish ice 
from liquid water on surfaces such as runways. In addition, the FLI may be able to detect 
icing conditions ahead of an aircraft during flight, by discriminating supercooled liquid water 
drops from ice crystals. 

The sensitivity and characterization studies reported here addressed the phenomenology 
that supports detection by an FLI. Techniques for determining the range, and hence warn¬ 
ing time, were demonstrated for several of the hazards, and a table of research instrument 
parameters was developed for investigating all of the hazards listed above. This work sup¬ 
ports the feasibility of detecting multiple hazards with an FLI airborne sensor, providing 
enhancements to aviation safety. Future research should also address the additional infor¬ 
mation available in two dimensional images. This could be of particular value for turbulence 
and wake vortex detection, and is essential for enhanced imaging applications. Once feasi¬ 
bility of detection and determination of hazard levels has been established, more attention 
should be focused on issues including probability of detection, false alarm rates, efficient 
real-time processing algorithms, generation of displays, and instrument considerations such 
as hardware specifications required for selected hazards and sensitivity to vibrations. 

FLI technology shows promise for a new class of airborne sensors; however, further re¬ 
search must be performed to develop a means to estimate the intensities of the hazards 
posed to an aircraft and to develop robust algorithms to relate hazard levels to sensor read¬ 
ings. In addition, validation tests need to be performed with a prototype system. In order 
to accomplish these goals, additional research should include ground-based measurements, 
airborne measurements, further sensitivity studies, and simulations. Additional possible 
applications not addressed in this study should be investigated. These include detection 
of runway contamination, observations of temperature and water vapor for use in weather 
models, and detection of pollution or other aerosols. Downward-looking atmospheric pro¬ 
filing for observations of quantities such as water vapor, temperature, and turbulence also 
merits consideration. 
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Appendix A 

Workshop and Results 


The Airborne Radiometric Detection of Aviation Hazards Workshop was held 14-15 
August 2006 at the National Institute of Aerospace (NIA) in Hampton, VA. The goals of 
the workshop were to review the current status of research in this technical area, to assess 
the resources that are available for further work, and to familiarize interested participants 
with the research program. The workshop had 30 attendees from government, industry, 
and academia, with 11 technical presentations and additional discussions. The Workshop 
agenda and presentations are available online at a website hosted by NIA: 

http://www.nianet.org/workshops/ardah06.php. 

The Workshop resulted in lists of resources and other parties interested in participating 
in future research efforts. Simulation resources are available at Hampton University, the 
University of Wisconsin-Madison, and NCAR. Measurement resources will be provided by 
industry. Possibilities for equipment include Telops, Bomem, ITT, and Dr. Fred Prata’s 
multi-spectral imager. Additional parties who will be included in the research as it moves 
forward are airline representatives, including Delta and Northwest Airlines, pilots, the Fed¬ 
eral Aviation Administration (FAA), and airframe manufacturers, such as Boeing. Hampton 
University’s modeling capabilities lie in single line-of-sight radiative transfer. Another po¬ 
tential resource is the University of Wisconsin-Madison SSEC/CIMSS (Space Science and 
Engineering Center/Cooperative Institute for Meteorological Satellite Studies) group, which 
has extensive modeling capabilities that include a 4-D simulation and visualization capabil¬ 
ity that can be used to generate modeled images of water and ice clouds and volcanic ash. 
ITT also has modeling capabilities that include simulation and visualization. 

NCAR is the leader in turbulence research and forecasting. They have a broad database 
of information that can be used in future modeling efforts. Some data has already been 
furnished to HU by NCAR for the purposes of this reported feasibility study. 

Telops has a lightweight and compact hyperspectral imaging sensor using FT-IR tech¬ 
nology called FIRST. Other sensor technology options include Prata’s multispectral imager 
that is currently configured for volcanic ash detection, ITT’s interferometric imager, and 
Bomem’s lamellar grating Fourier Transform Spectrometer system. Prata’s imager is a 
non-commercial product that is configured for volcanic ash and SO2 detection. It has been 
tested on the ground at several volcanoes. Each system has its own characteristics, and 
more than one system may be warranted for ground/flight tests, keeping in mind physical 
and budgetary limitations. 

Other interested parties include airline representative from Delta and Northwest. North¬ 
west Airlines was unable to attend the Workshop, but would like to participate in the 
research as it moves forward. Delta Airlines stated that their interest in CAT is mainly 
economic. Delta, as an example of a commercial airline that would be interested in mount¬ 
ing the sensor that results from the study into their fleet, understands the relevant needs 
and the business case for such a system. Lastly, Delta may offer a low-cost platform for 
evaluation of a commercial sensor. 
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Appendix B 


Additional Information on Low Level Wind Shear 
studies by Adamson 


Adamson has described airborne IR hazard sensor systems in three US patents [34], 
[35], [31]. The first device was a mechanically chopped single-band radiometer operating 
in the 20-40 /jm region. The detector was not cooled. Automatic calibration was achieved 
by periodically inserting a paddle heated to 42 °C into the path of the received radiation. 
Much of the patent is concerned with an analog signal processing scheme to enhance the 
signal-to-noise-ratio, but the frequency response and noise figures are not given. An earlier 
patent by Peter Kuhn is included by reference [36]. 

Adamson’s second device includes multiple bandpass filters on a rotating wheel. The 
filters span the wavelength ranges 14 -14.5, 14.5 -15, 15 -15.5, and 15.5 -16 /urn. The input 
radiation was mechanically chopped and the detector was not cooled. The signal processing 
electronics had a low-pass filter with a cutoff frequency of about 1.6 Hz. No noise figures 
are given. Again, much of the patent is concerned with an analog signal processing scheme 
and the patent by Peter Kuhn is again included by reference. 

Adamson’s third patent describes a sensor based on an IR grating spectrometer, again 
with built-in calibration. The spectrometer provides a bandpass with adjustable width and 
center wavelength. A scheme is presented for integrating the radiometer with other sensors 
in order to combine predictive (radiometer) data with reactive (accelerometer) data, and 
the severity of the hazard is quantified in terms of the F factor. 

Adamson named his sensor the Advance Warning Airborne Sensor (AWAS) and de¬ 
veloped three generations of it, known as AWAS I, II, and III. The AWAS sensors pointed 
straight ahead in the flight direction; they did not scan in angle. During the summer of 1990, 
AWAS III was installed on a Cessna Citation operated by the University of North Dakota 
and operated on 66 test flights that included wet microbursts in Orlando, Florida and dry 
microbursts in Denver, Colorado. Good correspondence was obtained between hazard values 
predicted by AWAS and data from a Terminal Doppler Weather Radar (TDWR). 

During February - June 1991, TPS installed AWAS III sensors with data recorders on 
three MD-80 aircraft operated by American Airlines, and three more on DC-9 aircraft 
operated by Northwest Airlines. The sensors recorded data during more than 10,000 hours 
of commercial flights. One important finding was that the system generated a “nuisance 
alert” on 20% of the American Airlines flights. 

The NASA flight tests were reported by McKissick at the 5th wind shear conference [37]. 
With funding from a NASA Phase I SBIR grant in 1987, TPS had determined that a passive 
IR system appeared to be feasible for wind shear detection. Under a Phase II SBIR (1989- 
1991), AWAS I was flight tested on the NASA 515, and AWAS III was developed. 

AWAS III was then flight tested on the NASA 515 at Orlando and Denver in 1991 and 
numerous changes were made, including improvements in the optics and periscope, a new 
method for compensating lapse rate effects, and filtering of the hazard indices. The hazard 
indices were applicable for approach speeds only (approximately 140 knots) and NASA’s 
flight tests were conducted at much higher speeds, so an approximation was developed 
by NASA personnel to transform measured F-factors for comparisons with the AWAS III 
hazard indices. The improved AWAS III was flight tested in 1992, again on the NASA 515 
at Orlando and Denver. 

NASA’s conclusions from the AWAS III flight tests were that the sensor did provide 
a predictive capability for wind shear with 11 to 55 seconds warning time, but that the 
hazard indices were not reliable predictors of measured F-factors, and so the development 
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of a reliable hazard index was incomplete. At the 4th wind shear conference, Adamson 
discussed the reasons why the NASA flight test results were not as promising as the previous 
results [38]. He stated that the NASA test had the following problems: 

1) NASA personnel manually overrode the AWAS automatic mode selection routine 
(TPS personnel were not on the flights). 2) The location and method of mounting the IR 
window assembly impaired the sensor’s line of sight and led to accelerated mirror damage 
that was not seen on previous flights, as well as problems with rain. 3) Airspeeds during 
the tests were in excess of the normal landing and takeoff speeds that the AWAS algorithms 
were developed for, and NASA’s approximations for adjusting the test data to compensate 
for this problem were inadequate. 

In a question-and-answer session after Adamson’s talk, the consensus of the NASA par¬ 
ticipants appeared to be that IR radiometric technology was not sufficiently robust to provide 
both acceptable detection and false alarm rate performance. 

NASA’s previous conclusions about radiometry pertain to the AWAS series of sensors, 
which were all single-line-of-sight radiometers pointed along the direction of flight. The 
temperature resolution of the AWAS radiometers was not given in the patents nor in the 
conference proceedings, but it is perhaps on the order of ± 1 °C due to the use of un-cooled 
detectors. The spectral resolution of AWAS was not given either, but the radiometer with 
multiple filters on a wheel used 0.5 /mi bandwidths. 

In contrast, the FLI investigated in the study reported here is an imaging system that 
will provide information on the structure of the disturbed air ahead of the aircraft, with 
a temperature resolution on the order of 0.1 C and spectral resolution on the order of 0.5 
cm -1 (0.005 /tm at 10 /im wavelength), which is sufficient to resolve individual lines in the 
CO2 spectrum and the windows between them. The combination of high spatial, spectral, 
and temperature resolution will provide a wealth of detailed data that has never before been 
obtainable, and that data is expected to enable much more sophisticated algorithms with 
high detection rates and low false alarm rates. 
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Appendix C 


Further Discussion of Candidate Sensor Technologies 


Discussions have been held with the Telops Corporation in Quebec, Canada and ITT 
Industries Space Systems Division in Fort Wayne, IN. Telops produces the FIRST instru¬ 
ment, one of the candidate imaging interferometers for experimentally validating the FLI 
aviation hazard phenomena by the detection techniques outlined in the report. ITT has 
a prototype hyperspectral testbed system, originally developed for the Hyperspectral En¬ 
vironmental Suite (HES). Each of the companies’ instruments is discussed in more detail 
below. 

C.l Telops 

Telops’ FIRST sensor is shown in figure Cl [39]. FIRST is a lightweight and compact hyper¬ 
spectral imaging sensor using FT-IR technology. It provides unparalleled spatial and spec¬ 
tral information about the IR targets under measurement. With its integrated calibration 
system and electronics, FIRST can be used in a standalone (and unmanned) configuration 
for over 24 hrs. The standard hardware interfaces make the sensor easy to use. 



Figure Cl. Telops’ FIRST Imager. 

The sensor specifications are nearly ideal for the purposes of the planned data collection 
efforts if the wavelength range can be expanded to include 12-15 microns. The system is 
compact and portable for field use and has a wide field-of-view and high spatial resolution, a 
variable spectral resolution down to 0.25cm _1 , flexible data products utilizing real-time on¬ 
board digital processing, integrated calibration targets that ensure accurate measurements, 
simultaneous visible imaging, and an easy-to-use FT Pro data management suite. Table Cl 
fists the specifications of the FIRST imager [39]. 

However, the system operates in the 7.8 — 11.6pm range OR the 2 — 5.5pm range. One 
system cannot operate in both bands. Telops owns a long-wave unit for demonstration 
purposes only. Telops does not have a mid-wave unit in hand, but they have two customers 
who own mid-wave systems. These customers are the University of Madrid and Air Force 
Institute of Technology (AFIT) in Dayton, OH. 

Making modifications to the system is not trivial. Extending the long-wave range below 
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Table Cl. TELOPS’ FIRST IMAGER SPECIFICATIONS 


Unit 

Specification 

Spectral Band 

lim 

8-11 (LW) 3-5 (MW) 

Number of Pixels 


320 x 256 

Spectral Resolution 

cm -1 

0.25 - 150 

IFOV 

mrad 

0.35 

Communication 


Ethernet 

Data Transfer 


Cameralink 

Acquisition Software 


FT Pro 

Detector Cooling 


Closed Cycle 

Power Source 


28 V (DC) or 129 V (AC) 

Weight 

lbs 

< 60 


7.8/j.m means changing a cold filter (this modification implies a custom instrument - there 
is no access to the cold filter). Perhaps the system can be extended to 6.5 /im, but the data 
will be noisier. To be beneficial to a follow-on project, the resolution would have to be on 
the order of a wavenumber. FIRST can achieve 1/4 wavenumber or higher resolution. The 
applications also require at least one frame per second (fps). The long-wave FIRST can 
achieve one fps, but only with a reduced image size. The mid-wave version has a higher 
frame rate [40]. 

Telops has conducted several tests and data collection efforts with their system. En¬ 
hanced visibility work has not been performed with this system, nor have flight tests. 

C.2 ITT 

ITT has developed a lab prototype hyperspectral system with internal research and devel¬ 
opment (IRAD) funding. The prototype was developed as part of a hyperspectral testbed 
program for Hyperspectral Environmental Suite (HES). HES has been canceled, freeing up 
the laboratory prototype for other uses, which potentially includes detection of aviation 
hazards 01 . 

Currently at ITT, further IRAD funding is being sought for calendar year 2007 to package 
the lab prototype system. Hampton University has provided specifications to ITT for the 
proposed aviation hazard detection system to support the ITT IRAD proposal. By fall of 
2007 through IRAD funding, the lab prototype may be turned into a system that could be 
used in the field and that meets the required specifications for the aviation hazard detection 
project. The current system operates from 3.5 to 16 /jm with a 0.5 to 1 wavenumber 
resolution. The whole imaging array is populated, and the data system is a computer 
cluster. The field of regard is somewhat small, and a telescope is needed. The telescope 
addition would also be funded under the IRAD. Currently, the system is three separate 
arrays. Each array has been procured and demonstrated individually. The system was 
designed to operate all three at the same time, and the ITT IRAD funding would provide 
the means to do the integration [41]. 


cl Glen Davis at ITT led the development effort for the HES testbed. 
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Airborne Forward-Looking Interferometer Turbulence 
Investigation 


Larry B. Cornman 

National Center for Atmospheric Research 
Boulder, CO. 


D.l Introduction 

Turbulence is the leading cause of non-fatal accidents for Parts 121, 135, and 91 air traffic. 
In addition, the costs to the major US commercial air carriers are well over $100 million 
dollars per year. NASA and the FAA have sponsored very successful programs to detect 
and forecast turbulence. NASA has focused on forward-looking airborne sensors, such as 
Doppler radars and lidars - and is currently is supporting a program to investigate the 
feasibility of using a passive infrared sensor to detect turbulence ahead of an aircraft. 

This report is divided into two parts: Part A discusses the work done investigating 
methodologies for measuring turbulence from an airborne IR interferometer; Part B provides 
an overview of the data sets that were delivered to Hampton University in support of their 
simulation activities. 

Part A 


D.2 Theory 

It should be noted that the material presented below reflects preliminary work. A number 
of assumptions are made, which need to be supported. Furthermore, the equations need to 
be checked for errors, as well as further simplifications where possible. Nevertheless, it does 
present the basic concepts and methodologies. 

Following on the methods used in the airborne radar problem, the concept is to derive 
equations that relate the statistics of the atmospheric turbulence (e.g., temperature field) to 
those of the sensor measurements (e.g., the irradiances). Two approaches can be considered: 
(1) the turbulent temperature field is homogeneous and extends along the whole measure¬ 
ment path, and (2) the temperature field is laminar, excepting for a ’’patch” of turbulence 
at a certain distance from the sensor. 

Start with the Radiative Transfer Equation for a gas in thermal equilibrium, with neg¬ 
ligible scattering, and the measurement is made along the direction x: 

dI ^'= h(x)\W B (x) - I[x)\ (Dl) 

where I(x) is the irradiance (or intensity), Wb(x) is the Planck function and h{x) is the ab¬ 
sorption coefficient. It should be noted that the frequency dependence of each of these terms 
is being suppressed for notational clarity. Equation (Dl) is a linear first order differential 
equation of the form, 


dy _ 

dx 


a(x)y(x) = b(x) 


With the solution, 

y(x) = exp J o(a;)da;| J b(x) exp J a(a; , )dx , | dx + const. 


(D2) 


(D3) 
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Applying this to equation (Dl), gives: 


I(L) — 1(0) = exp J ft(a;)dx| J H(x)Wb(x) exp J h(x , )da; , | dx (D4) 

e., there is no radiation coming fro 

/ h(x')dx' = — h(x')dx‘ 

Jx Jo Jo 


Assume that I(L) is zero, i.e., there is no radiation coming from x > L. Further, note that 


(D5) 

and since the integral over [0,L] is independent of x, it can be taken outside the integral 
over x, which gives 


1 (0) = — J 1i(x)Wb(x) exp J /i(x , )dx , | dx (D6) 

This is a well-known result in the literature (e.g., Tourin and Krakow, 1965), the full deriva¬ 
tion was given above to fill in the steps. Let (3(x) = / x k(x')dx\ the so-called optical depth. 
Note that 


and hence 


d(3(x) 

dx 


= h(x) 


(D7) 


i exp [-1 ^ exp i-'** 1 *) <D8) 

Hence, equation (D6) can be written as 

1(0) = J W B (x) exp [-p(x)\ dxi^j W B (x) exp [-/3(x)]d/3(x) (D9) 

This is a well-known form given in the literature (e.g., Gautier and Revah, 1975). Another 
way to write equation (D6) is to let a(x ) = exp [—(3(x)], the transmittance. Equation (D9) 
then becomes 


7(0)= [ W B (x)-^a(x)dx = [ W B (x)da(x) (DIO) 

Jx =0 dx J x — o 

It should be noted that some authors use the notation t(x) = /3(x), while others use t(x) = 
a(x), so care should be taken when trying to compare references. 

Consider another measurement made at the point, x = p > 0, 




h(x)WB(x) exp 


[-/; 


h(x r )dx r dx 


This means the correlation between 7(0) and I(p) is 


(Dll) 


/ pL pL+p px px' +p [ 

(I(0)I(p)) = C 2 (J J h(x)h(x')W B (x)W B (x')ewp - I -J h(y)dy\dx'dx 

(D12) 

There are references that derive results for the mean value of equation (Dl) for turbu¬ 
lent conditions. These include: Kabashnikov and G.I. Remit (1979), Kabashnikov (1985), 
and Kabashnikov and Myasnikova (1985). These references deal with correlations between 
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the absorption coefficient and the Planck function. These references discuss assumptions 
that can be made to simplify the fairly intractable equations. Two references have been 
found that deal with the correlation (/(O)J(p)), Rybicki (1965, 1967). They do not take 
into account correlations between the absorption coefficient and the Planck function, and 
only consider the fluctuating parts of these two functions. In the following, we will compute 
( I(0)I(p )) - assuming no correlation between the absorption coefficient and the Planck func¬ 
tion - but consider both the mean and fluctuating components. The increase in complexity 
in this approach is significant. Further work will have to deal with simplifying - yet justified 
- assumptions that can make the problem more tractable. Allusions to this process will be 
made below. 

Assuming that the Planck function satisfies the Rayleigh-Jeans Law (low frequency 
limit), Wb{x) = CT(x), where C is a constant (though a function of the frequency). De¬ 
compose the quantities into a mean and fluctuating part, T = (T) + V and h = (h) + h', 
which by definition means that (T) = 0 and (h) = 0. Using the fact that ((Z)) = (Z), the 
integrand of equation (D12) can then be computed. 


Where, 


and 


[T'(x)T'(x') + (T(x)) (T(x'))]e 


-Wl(X,X ,p) e ~LJ2(X,X , P ) 

(D13) 

+ ti(x) {h{x'))+ti{x)ti{x')\ 

\ [ + [ (h(y))dy\ 

(D14) 

[Jo Jp J 


\l" ■ f X +P (h'(y))dy\ 

(D15) 

[Jo Jp 



For the moment assume that the fluctuating parts of the absorption coefficient and the 
temperature are statistically uncorrelated. Taking the ensemble average of equation (D13) 
gives 


[(T'(x)T'(x')) + (T(x)) (T(x'))} ©! 

• [(h(x)) (h(x ')> (62) + (h(x)) (h'(x')Q 2 > + (h'(x)Q 2 {h(x ')» + {ti(x)h'(x')e 2 )] 1 J 

This is a rather complicated expression, with double- and triple-correlations of the h' func¬ 
tion. For now, consider the term, ^ e -Qi{x,x ,p)^ jf fo' j s assumed to be a Gaussian random 
variable, then 

^ e -<W»y,p>y = exp ( D17 ) 

where, 


={Qi{x,x', P ))=\r f +2rpp 

J 0 Jo Jo Jrho Jrho 


{h'{y)h'{y'))dydy'\ (D18) 


Note that for a homogeneous h' field, (h r (y)h'(i/)) = R h >(y' — y). In the first double integral 
make the change of variables y = y' — y, y = (y 1 + y) /2. This gives 


a: 


Rw ( rj)drjdy 


(D19) 
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If it is assumed that Rh(v) —> 0 for rj » x, then the integral over pin equation (D19) can 
be written as 


f Rh'(v)dv « f Rh'(v)dr] = 2 L w (h ,2 (p)) (D20) 

Where the integral scale of h', L h >, is defined by 

Lh ’ = WW)l Rh ' ( ' l)d ' 1 (D21) 

And so, equation (D19) becomes 

f X [ X Rh'(rj)di}dii * 2 L h , f ( h ,2 (p)> d,u (D22) 

Jo J-X Jo 

In a similar fashion, the other two double integrals in equation (D18), are given by 


nX rx'+p rX+X+p /2 

2/ / {ti{y)ti{y')) dydy' = 2L h , / (h' 2 (p))dp (D23) 

JO Jp J p/2 


rx'+p rx'+p rx+x'+p /2 

/ / (ti(y)h'(y')) dydy' = 2L h > / ( ti 2 {p))dp (D24) 

Jp Jp Jp /2 

So finally equation (D17) becomes 


(' e-M*’*'* )) = exp |i h , |jT + 2 +P/2 + J* +P/2 | (D25) 

At least we can now evaluate the terms in equation (D16) - excepting for the double- and 
triple-correlations of the function (that will have to come at a later date). Recall that these 
terms are the integrand of equation (D12). Let us call this (7(0)7(p)) 1 


f L fL+p ( px rx'+p 

{mmfi = C 2 J 0 J -x) + {T{X)) (T(x')}] expl -[jf +J (h(y))dy] 

■ (h(x)) (h{x')) exp\\L h ' +2 J + J jdxdx' 

(D26) 

Assuming homogeneous turbulence, Rt{x' — x) = C'^gi/x' — x) ■ Cf. is the intensity pa¬ 
rameter of the turbulent temperature field - also known as the structure constant of the 
temperature turbulence, and g(x' — x) is an assumed known functional form (e.g., Kol¬ 
mogorov or von Kdrmdn). Kabashnikov, (1985), in studying the equation for (7(0)), dis¬ 
cusses conditions (the so-called optically thin fluctuation approximation, or OTFA) under 
which Then equation (D26) is 


pL nL-\-p 

(7(0)7(p)) 1 p C 2 C 2 , J J [g(x' -x) + (T(x)) (T(x'))} 
■ exp/— [/ + J (/i(p))dp| 1 (h(x)) (h(x')) dxdx' 
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In principle then, if the mean values of h can be estimated (from, say, a standard atmosphere 
model), and if the h-correlation terms in equation (D16) are small relative to the other 
terms (remains to be shown), then from equation (D27) the structure constant Cj., can be 
estimated by the ratio: 


_ mm _ 

° 2 Jo Jp +P - x ) + ( T ( x '))} exp{- [/ 0 X + Jp +P {h(y))dy\ } ■ (h(x)) (h( x')) dxdxJ 

(D28) 

The denominator term can be computed a priori as a function of altitude, with the as¬ 
sumption that the airplane is in straight and level flight. (J(0)/(p)) would be calculated 
from the irradiance measurements at a given frequency. Evoking the Ergodic hypothesis, 
(replacing ensemble averages for space averages, given a homogeneous turbulent field), the 
average value, ( I(0)I(p )), would be computed as a running spatial average (i.e., as the air¬ 
craft moves). Furthermore, instrumental characteristics (as a function of frequency) can be 
accommodated. Finally, since equation (D28) is an estimate for the structure constant at 
a given frequency, an averaged value of the structure constant can be constructed via an 
average over frequency. Note that the structure constant is frequency-independent, all the 
frequency dependence is on the right had side of Eq. (28). 

Another approach - using the same methodology as presented above - would entail the 
development of a similar relationship between the power spectrum of the irradiance measure¬ 
ments and the spectrum of the temperature field. This would be accomplished by taking the 
Fourier transform of equation (D27). A third approach would be to use structure functions. 
The second-order structure function for the irradiance field is given by 

Di{ 0, p) = ([I{p) - /(0)] 2 } = </ 2 (p)> + </ 2 (0)> - 2 </(0)/(p)> (D29) 

We note that the last term in this equation is given directly by equation (D27). The first 
two terms of equation (D29) are given by replacing I (0) by I ( p ) and vice-versa, respectively. 
That is 


and 


(7 2 (p)) 1 s C 2 C 2 , j J W - x) + (T(x)) (T(x'))} 
■ exp/ — [/ + J (/i(y))dj/| 1 (h(x)) (/i(x')} dxdx' 


(D30) 


(I 2 (0)) 1 w C 2 C 2 £ £ [g{x’ -x)+ (T(x)) 

■ exp/— [ J + J (/i(p))dy| 1 (h(x)) (h(x')) dxdx 1 


(D31) 


It is not clear whether using the structure function approach will result in any cancellations 
of the terms in equations (D16), (D27), (D30), and (D31). Simplifications often occur in 
the application of structure functions, especially a separation between mean and fluctuating 
parts, however, further analysis is required to flesh this out. 

The methodology presented above assumed that the turbulence was homogeneous over 
the path (0,L). However, it can also cover the case of inhomogeneous turbulence with slowly 
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varying intensity (e.g. structure constant). As mentioned above, another avenue to investi¬ 
gate is the model of a fairly laminar atmosphere with a high-intensity, small patch (1-3 km 
wide, say) of turbulence a certain distance from the aircraft. It should be noted that this 
work is in a very preliminary stage; however, a sketch of the concept is given below. 

Following Tourin and Krakow (1965), though reversing directions in their derivation, 
we consider taking the atmosphere along the sensor pointing direction and breaking it into 
n-distinct regions. Then 7(0) = Y^i=i where /,• is the x-value for the i th region. From 
equation (D6) 


1(h) = - j' h(li)W B (li) exp fe(a; / )da; / j exp — ^ j h(x')dx' 


Working through all the steps, and letting «(//,) = exp—h(lk)[lk — h-i], yields 


(D32) 


So, I (0) is given by 


m = w B{h) i<m i] 


(D33) 


7(0) = J2 W b(U) Hh) - 1] n (D34) 

Now assume that all the regions i = 1,7 — 1 and i = k +1, n are laminar and i = k is 
turbulent. As before, let W B (l k ) = CT(l k ), and T = (T) + T' and h = (h) + h!, this gives 


I(l k ) = l(T(l k )) + T'(l k )] M4) -1] n mlh) (D35) 

j=o 

where 

ai (l m ) = exp{- [(h(l m )) + ti(l m )\ [l m - l m - 1 ]} (D36) 

Now the rest of the summation in equation (D34), will just contain the average quantities 
which can be estimated using a standard atmosphere (same for the average quantities in 
equation (D35)). The next steps could follow a similar path as that shown above: compute 
correlation, spectral or structure function equations - or even just computing the mean- 
square value of 7(0). These paths will have to be examined during any follow-on activity. 


D.3 Estimating the impact on an aircraft 

Unfortunately, the atmospheric quantity that dominates the response of an aircraft to tur¬ 
bulence is the vertical component of the wind vector - not the scalar temperature field. 
Therefore, a theoretical or empirical relationship between these two quantities must be 
established; specifically the connection between the intensity metrics, Cf. and the eddy dis¬ 
sipation rate to the two-thirds power, e 2 / 3 . Given a von Karman form for the velocity field, 
a straightforward relationship between e 2 / 3 and the variance of the velocity field, cr 2 ,, can be 
made. Furthermore, in the airborne radar program, the relationship between the standard 
deviation of the velocity field (the square root of the variance) and the root-mean-square 
(RMS) vertical accelerations of a given aircraft was developed. Therefore, in principle, the 
IR sensor measurements can be related to the impact on an aircraft. 



Of course, the term ”in principle” hides the fact that a lot of work is required to show 
that these concepts are applicable to the real world and that a system can be developed 
(sensor and algorithms) to produce operationally useful performance metrics, e.g., probabil¬ 
ities of detection (PODs) and false alarm ratios (FARs). In the airborne radar program, it 
was found that analytical and realistic simulations of the turbulent field and the subsequent 
sensor measurements were invaluable in the development and testing of the turbulence de¬ 
tection algorithms. Of course, simulations can go only so far in this process; hence, field 
testing is another important component of a successful program. 


D.4 Feature extraction from the sensor data 

Another aspect of detecting turbulence from a passive IR sensor is feature extraction. This 
is in distinction from the turbulence intensity methodology described above. By ’’feature” 
we mean characteristics in the sensor data that are the result of the turbulence - and not 
the background. Preliminary feature extraction algorithm concepts have been applied to 
simulated brightness temperatures and interferogram data. 


D.4.1 Brightness temperature data 

Consider the severe turbulence encounter case shown in figure D1 below. The top panel in 
this figure is the vertical acceleration time series and the bottom panel is the temperature 
field. Figure D2 illustrates the brightness temperatures differences calculated from the IR 
sensor simulation. That is, the measured temperature field from the encounter is used 
to simulate the irradiances, which are then converted into brightness temperatures. The 
brightness temperature is computed as the temperature that satisfies, T n = IT/?(/((])). 
The brightness temperature difference is obtained by subtracting the standard atmosphere 
temperature from the brightness temperature. 
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Figure Dl. Vertical acceleration and temperature time series. 
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In figure D2 the black curve is the brightness temperature differences. The red curve is 
obtained by using a wavelet filter. Wavelet transforms decompose the signal into coefficients 
that are a function of two parameters: the “scale” parameter and the “position” parameter. 
As a comparison, the Fourier transform decomposes the signal into coefficients that are a 
function of the single parameter, wavelength. The scale parameter in the wavelet transform 
is akin to the wavelength parameter in the Fourier transform. The position parameter, in this 
case, is the wavenumber value. The advantage of wavelet transforms is that discontinuities 
in the data can be resolved without “smearing,” as would occur with a Fourier transform. 
The wavelet filter used to produce the red curve in this figure is obtained by setting the 
coefficients above a chosen scale value to zero, and then reconstructing the signal. It can be 
seen that this wavelet filter acts in part like a Fourier low-pass filter, removing the higher 
frequency oscillations, but unlike a Fourier filter it leaves the larger-scale discontinuities 
intact. 

The wavelet filter is useful in helping to isolate the discontinuities from the higher-scale 
regular oscillations; however, it does not determine which points are associated with the 
discontinuities. A valuable tool for performing this task is the UCAR-developed Intelligent 
Outlier Detection Algorithm (IODA). Referring to figure D2, it can be seen that the dis¬ 
continuities in the wavelet-filtered signal can be thought of as outliers or distinct features 
separate, from the rest of the wavelet filtered signal. IODA is a powerful algorithm for fea¬ 
ture detection, and as seen in figure D3 - applied to the wavelet-filtered signal - it isolates the 
’’regular” points (inside the contours) from the discontinuities (points outside the contours). 
Thus, figure D3 is the same data as shown in figure D2 but on a different scale, illustrating 
how IODA identifies those points associated with the sharp discontinuities (points outside 
the contours). 


D.4.2 Interferogram data 

Another prototype feature-extraction algorithm was created for the interferometric data. 
This process contained a number of steps. The first step was to take the raw interferogram 
and compute local minima and maxima, in effect creating an envelope function. This can 
be seen in the top panel of figure D4, where the black points are the raw data and the red 
points are the local minima and maxima. It is easy to see that the envelope of the data 
between optical depths (OPD) 0.6 and 0.7 cm is quite different that that of the rest of the 
data. That is, the envelope in this region is very monotone, as opposed to oscillatory in 
the other regions. This naturally leads to the concept of creating a ’’difference” function, 
the difference being the local maxima subtracted by the local minima. This is shown in the 
bottom panel of figure D4. There is clearly an exponential trend in the difference function, 
so the next step was to remove this trend. This was accomplished by doing a line-fit in 
log-space, between approximately OPDs 0.05 and 0.55 - as can be seen by the red line in 
the second panel from the top of figure D5. The trend curve is shown in linear-space in the 
third panel from the top of figure D5. Finally, the bottom panel of figure D5 shows the 
effect of the trend removal. Since this trend-removed data is still somewhat noisy, a wavelet 
filter was applied to it, and the results are shown in figure D6. The oscillatory nature of 
the data at OPDs out to approximately 0.6 cm is quite clearly seen. Given the available 
time, this was as far as the methodology was worked out. However, it is clear that the next 
steps involve the differentiation between the oscillatory parts and the non-oscillatory parts. 
A preliminary attempt, using running mean filters, was made - but with mixed results. 
Nevertheless, we are confident that the final feature detection can be made. Obviously, 
more cases need to be examined, especially cases where the turbulence was further from the 
aircraft. At those distances, the signal becomes weaker, so having an algorithm that can 
extract features associated with the turbulence becomes even more important. 
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Figure D3. IODA identification of sharp discontinuities. 
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Figure D4. Interferogram, envelope function, and envelope differences. 
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Figure D5. Difference function, trend line, trend curve, and difference signal. 
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Figure D6. Wavelet filter applied to trend-removed difference function. 
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Part B 


D.6 Inputs to FLI simulation 

The FLI produces both the interferogram and the retrieved brightness temperature from 
the temperature and humidity profiles along the viewing axis. For airborne detection of 
turbulence, the FLI would observe atmospheric conditions ahead of the aircraft with a fixed 
line of sight with a few different elevation angles. Accurate warnings of a turbulence hazard 
requires a reliable connection between thee FLI signatures and the turbulence in the vertical 
velocity field on the scales of motion that affect a given aircraft, typically 100-1000 m. In 
addition, the number of false alarms must be small to gain the confidence of the pilots. So 
there are two fundamental research tasks that must be performed as part of the development 
of a reliable FLI turbulence detection algorithm: 

• Develop methods to differentiate temperature signatures that are turbulent from non- 
turbulent spatial temperature variations, such as fronts and even waves, 

• Develop a (probably statistical) connection between temperature turbulence to veloc¬ 
ity turbulence. 

These tasks can be readily addressed in a simulation environment where realistic turbu¬ 
lent and non-turbulent atmospheric conditions are contained in a digital data cube which 
can be used as input to an FLI simulation. In initial work preformed in FY06 NCAR pro¬ 
vided Hampton University some turbulence data to be used as input to the simulations. 
These were in the form of three-dimensional idealized turbulence simulation volumes, and 
aircraft data from a commercial flight data recorder and a research aircraft, both of which 
experienced a strong turbulence encounter. 

D.7 Idealized turbulence simulations 

A first evaluation of the FLI system for turbulence detection was produced by extracting 
the FLI signatures from idealized high resolution numerical simulations of von Karman 
turbulence. In spite of some limitations, the von Karman turbulence spectrum has been 
shown to provide a reasonable approximation for localized turbulence events at mid to upper 
levels (Murrow et al. 1982, Murrow 1987). 

D.7.1 Background 

Assume that each component of the turbulent velocity field is described by a homogeneous 
and isotropic Gaussian random process. Then the turbulence statistics are completely de¬ 
fined by the longitudinal ( Bll ) and transverse (Bjvjv) covariance functions given by (Monin 
and Yaglom 1975, p. 19) 


Bll{t) = ^ [u(a;) - (u) ] [u{x + r) - {u) ] ^ 

(D37) 

B NN (r) = ^ [v(x) - {v) ] [u(x + r) - {v) ] ^ 

(D38) 


where u(x) and v(x) are the velocity components along and transverse to the displacement 
vector r = ( x,y,z ), respectively, and () denotes an ensemble average. The longitudinal 
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(Fl) and transverse (F/v) one-dimensional spatial spectra of the velocity field are given by 
(Monin and Yaglom 1975, p. 43) 


^=-J e~ ikr B LL (r)dr 

(D39) 

and 


Fjv(fc) = — J e~ lkr B NN (r)dr 

(D40) 

For the von Karman turbulence model (e.g., Hinze 1959) 


a 2 L 

~ 7T(1 + F 2 fc 2 ) 5 /6 

<7^(3+ 81^) 

FN{k) - 677(1 + L 2 fc 2 )iV6 

(D41) 

(D42) 


where a 1 is the variance of the given velocity component, L 0 is the outer scale of turbulence, 
and 


Li = = 0.7468F o (D43) 

is the integral length scale. The approximation for high wavenumber k (Monin and Yaglom 
1975, p. 355) is 

9 ( 0 , 2 / 31 . 5/3 

F L {k) = 3r(1/3 = 0.4977096c 2 / 3 A: -5 / 3 (D44) 

0(0 2/31.-5/3 

F N (k) = Qr(1/3 = 0.6636128e 2/3 A;“ 5/3 (D45) 

where e is the energy dissipation rate, and the Kolmogorov constant C is taken as 2.0 (Monin 
and Yaglom 1975 pp. 483-485). Equating the high wave number form of equation (D41) 
and equation (D44) produces 


L 0 = 0.933668cr 3 /e (D46) 

which connects the outer scale to the variance and energy dissipation rate. 

The spatial statistics of the temperature field for localized turbulent events can also 
be approximated by the von Karman model for the longitudinal velocity, i.e., the one¬ 
dimensional spatial spectrum is 


W'/idr < D47 > 

and the inertial range behavior becomes 

9(^2 u- 5/3 

F T (fc) = 3 g (1/3 = 0.2488548C 2 AT 5 / 3 (D48) 

where (7f. is the temperature structure constant. 

The spatial spectra F q (f) can be related to the corresponding temporal spectra S q (f), 
S v ( f ) and S w (f) derived from in situ tower data or aircraft data assuming Taylor’s frozen 
hypothesis is valid (Hill 1996; Wyngaard and Clifford 1977). The temporal spectra become 
(with k = 27 rf/U ) 
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s q (f) = 2 ij F <i( 2 l f ) ( D49 ) 

where U is the magnitude of the mean along stream velocity or the true airspeed of an 
aircraft. 

D.8 Simulation of von Karman turbulence 

One-dimensional homogeneous and isotropic turbulence fields can be generated for a wide 
range of turbulent conditions by constructing many realizations of a von Karman tem¬ 
perature field. Each realization is generated with a spectral technique that produces the 
exact spatial covariance defined by the parameters of a von Karman model (Frehlich 1997), 
i.e., the random temperature is produced by generating statistically independent zero mean 
Gaussian random numbers for the real and imaginary parts of the Fourier coefficients. The 
variance of each Fourier coefficient is chosen to produce the desired spatial correlation func¬ 
tion. The spatial realizations are calculated from the random Fourier coefficients using the 
Fast Fourier Transform (FFT). Each simulation produces two (one-dimensional) turbulence 
fields - one corresponding to the real part and the other to the imaginary part of the FFT. 
One simulation output (with to two realizations) was provided to Hampton University on 6 
June. Each realization contained ID temperature fluctuations in the longitudinal direction 
with a spacing of 1 m. This was produced using an outer scale of 2000 m and a standard 
deviation of 1 K, which gives a peak-to-peak variation of 4 — 5 K. The output could be 
scaled up or down to correspond to different turbulence intensities. These realizations of 
temperature were scaled to typical temperature signals for KH instabilities observed in re¬ 
search aircraft measurements (Whiteway et al. 2004) and located at various distances in 
front of the aircraft. 

The simulation algorithm is numerically efficient and can be easily extended to two- and 
three-dimensional random processes (Frehlich et al. 2001). In principle, the von Karman 
model for the temperature and velocity fields could be modified to include correlations 
between the temperature and velocity. However, this would require more research to provide 
a defendable correlation functional form. 

D.9 Research Aircraft Data 

Since a von Karman field does not provide a realistic connection between the temperature 
and velocity fields, examples of traces of spatial variability derived from the NCAR research 
aircraft were collected and provided to GTRI/HU for evaluation as ASCII data files. The 
test cases are described below. 

D.9.1 INDOEX 

The NCAR C-130 research aircraft deployed as part of the INDOEX (Indian Ocean Exper¬ 
iment) campaign in the tropics during 1999 contained a turbulent event (see figure D7) at 
an altitude of 4.8 km. This figure shows an INDOEX case with time series plots of along 
track velocity u, across track velocity v, vertical velocity w, temperature T, true airspeed 
(TAS), altitude, heading (HDG), and water vapor density p as a function of track distance. 
Data file is indoex.l. The magnitude of turbulence of interest to aircraft response is defined 
by the spatial spectrum of vertical velocity S w ( k ) as a function of the spatial wavenumber 
k which is shown in figure D8. This figure shows spatial spectra of u, v, w, and T from the 
middle half of the data of figure D7. The A: -5 / 3 regions are indicated by dashed vertical 
lines. The classical k ~ 5 / 3 spectrum is produced from about fc=0.00.3 - 0.02 m -1 with an 
EDR value of approximately 0.30 rn 2 / 3 s -1 , which would put the turbulence experience in 
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roughly the moderate intensity category. The temperature and water vapor density have 
some correlated features to the enhanced vertical velocity and also to some features that 
are not correlated. This may be related to the conditions in the tropics. 


D.9.2 T-REX cases 

The new NCAR HIAPER G5 aircraft collected considerable data over the Rocky Mountains 
and the Sierra-Nevada Mountains in March and April of 2006 in support of the Terrain 
Induced Rotor Experiment (T-REX). This data was collected at a higher altitude than the 
INDOEX data. Case 1 is an example of a topographically-induced gravity wave (or mountain 
wave) recorded by the HIAPER is shown in figure D9 with the corresponding spatial spectra 
in figure DIO. Time series plots from T-REX case 1 with along track velocity u, across track 
velocity v, vertical velocity w, temperature T, true airspeed (TAS), altitude, heading (HDG), 
and water vapor density p as a function of track distance are shown in figure D9. Shown 
in figure DIO are the corresponding T-REX spatial spectra of u, v, w, and T from the 
middle half of the data of figure D9. The k ~ 5 / 3 regions are indicated by dashed vertical 
lines. The temperature fluctuations are large and should produce a large signature in the 
FLI analysis. However, the vertical velocity spectrum in wavenumber region fc=0.0025-0.01 
m -1 has an EDR of about 0.012 m 2 / 3 s _1 for the vertical velocity which is a rather low 
turbulence intensity for aircraft induced turbulence. This is a common condition in the 
upper atmosphere where most of the energy is confined to larger scales, i.e., fc=0.0001-0.002 
m _1 as shown in figure D8. 

Another mountain wave case (Case 2) observed by the HIAPER is shown in figure 
Dll with time series of along track velocity u, across track velocity v, vertical velocity w, 
temperature T, true airspeed (TAS), altitude, heading (HDG), and water vapor density p 
as a function of track distance. Figure D12 is the sampled spatial spectra of u, v, w, and T 
from the middle half of the data of figure Dll. The k ~ 5 / 3 regions are indicated by dashed 
vertical lines. This case is for a lower altitude and with a humidity feature nearby. Again, 
this is a case with a large temperature variation but an EDR of only about 0.0085 m 2 ' /3 s _1 
from the vertical velocity spectrum in the spatial wavenumber region that affect aircraft. 


D.9.3 ATReC case 

During the Atlantic THORPEX Regional Campaign (ATReC) the NOAA G4 weather re¬ 
connaissance aircraft experienced severe clear-air turbulence on 6 Dec 2003 while above 
a region of active moist convection, to the east of a surface low pressure system off the 
north-east coast of the United States. This case has been analyzed in detail and numerical 
simulations of the event have been performed (Lane et al. 2005). The in-situ measurements 
from the G4 were taken at 1 Hz time resolution, and are shown in figure D13. Shown 
in the figure are in-situ data from the NOAA G4 aircraft between 15:20:35 and 21:43:55 
UTC 6 December 2003 are time series plots of longitudinal velocity (u), transverse velocity 
(v), vertical velocity (w), temperature (T), true air speed (TAS), altitude (Alt), aircraft 
heading (HDG), and relative humidity (RH). The severe turbulence was initially encoun¬ 
tered at approximately 20:06.This figure shows that shortly after 2000 UTC, the aircraft 
measured strong fluctuations in all three velocity components, especially vertical velocity, 
temperature, true air speed, and relative humidity. These turbulent fluctuations occurred at 
approximately 13.5 km altitude. Of note is the relative humidity which attains a maximum 
value of approximately 40%. Therefore, this turbulence did not occur within cloudy air. 
Nevertheless signals are obvious in the relative humidity field, but not as well-defined in the 
temperature field. 

The in-situ measurements of velocity and temperature taken during 500 seconds sur¬ 
rounding the turbulence encounter were analyzed using spectral decomposition. The spec- 




tral analysis shows that all of these four variables possess a well defined inertial range at 
scales of motion less than approximately 1800 meters, with the inertial range extending to 
larger scales for the horizontal velocity components and temperature. Using the method de¬ 
scribed in Sharman and Frehlich (2003), the power spectra levels are used to determine the 
eddy-dissipation rate e 1 / 3 for the three velocity components and were determined to be e!/ 3 
= 0.301 m 2 / 3 s -1 (longitudinal), ei^ 3 = 0.277 rn 2 / 3 s _1 (transverse), and e^ 3 =0.139 m 2 / 3 s _1 
(vertical). In a similar way, Cf. was derived to be 0.029 K 2 m ~ 2 / 3 . The standard deviations 
of the velocity components and temperature are cr u =10.8 m/s, a„=8.0 m/s, a lu =().93 m/s, 
and ctt= 1.8 K, respectively. Figure D14 is a plot of longitudinal (u), transverse (v), vertical 
velocity (w), and temperature (T) spatial spectra developed from the time series of data in 
figure D13. The area contained between the vertical dashed lines exhibit a k ~ 5 / 3 spectral 
slope. Note that, in this case, the turbulent velocity fluctuations were much stronger in the 
horizontal components compared to the vertical component, emphasizing the anisotropic 
behavior of discrete events. 


D.10 Commercial aircraft Flight Data Recorder (FDR) 

Analysis of the FDR information from 7 turbulent events produced a couple of interesting 
cases. Unfortunately, in all cases, humidity information is not recorded, and for several of 
the cases temperature information was not recorded either. The three datasets that seem 
to have useful information are summarized below. 


D.10.1 Severe turbulence encounter case, Houston, TX, 6 Aug 2003 

This turbulence encounter case occurred on August 6, 2003 at 20:57 UTC as an Airbus 
A340 was at cruise altitude over Walnut Ridge, AR (approximately 36:33 N latitude, 90:42 
W longitude, and 31,000 ft (9448.8m)) en route to Houston, TX. The radar reflectivity (as 
measured by a nearby NEXRAD radar) grew from about 10 to 30 dBZ at the location of 
the encounter during this time, but its magnitude would likely not have appeared dangerous 
to the pilots. Figure D1 above shows the vertical acceleration and temperature time series 
from the flight data recorder. 


D.10.2 Convective case over southeastern Iowa, 22 Oct 2004 

A commercial aircraft encountered what was reported as severe turbulence over southeastern 
Iowa on 22 Oct. 2004 at FL370 at about 1636 UTC. A time series of the recorded vertical 
acceleration, total air temperature (TAT) and derived static air temperature (SAT) is shown 
in figure D15. This figure shows the FDR trace of turbulence encounter in or near convection 
over southeastern Iowa on 22 Oct 2004. Shown are the measured vertical acceleration (red) 
recorded at 8 Hz, TAT (green) recorded at 1 Hz, and corresponding SAT (blue) based on 
a constant Mach number during the encounter. The vertical acceleration is measured at 
8 Hz on the recorder and shows a large swing from a maximum of 1.80 g to a minimum 
of 0.18 g in just a few seconds. The TAT is recorded at 1 Hz and shows a sharp increase 
at about the time of the event, followed by a more gradual decrease and then followed 
by a gradual recovery. But the TAT includes the effects of adiabatic compression in the 
temperature probe. It is related to the SAT or ambient temperature through the relation 
SAT = TAT/( 1 + CM 2 ), where M is the aircraft Mach number and C is a sensor specific 
constant. For the cruise Mach number, SAT/TAT ns 0.5. Using this conversion gives the 
derived SAT curve in the figure. 
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D.10.3 Convective case over northwestern TN., 22 July 2002 

A commercial B757 aircraft encountered what was reported as severe turbulence over north¬ 
western Tennessee on 22 July 2002 at 1848 UTC at about 35,000 ft. while climbing to a 
destination cruise altitude of FL390. A time series of the recorded vertical acceleration and 
total air temperature (TAT) is shown in figure D16. This figure shows the FDR trace of tur¬ 
bulence encounter in or near convection over northwestern Tennessee on 22 Jul 2002. Shown 
are the measured vertical acceleration (red) recorded at 8 Hz and TAT (green) recorded at 
1 Hz. Abscissa is 8 Hz counts from 18:45:07 UTC. The vertical acceleration is recorded 
at 8 Hz and shows a large swing from a maximum of 1.74 g to a minimum of -0.37 g in 
just a few seconds. The TAT is recorded at 1 Hz. Unfortunately, SAT was not recorded, 
so a conversion factor could not be applied to derive a high rate SAT during the event. 
However, the TAT does show a gradual decrease, as expected during a climb, prior to the 
event, with oscillations coincident with the event. Thus in this case there seems to be a 
temperature signature during the time of the event, but it is difficult to know how much of 
this may be due to Mach variations. Note though that this event occurs only over a very 
small time/space interval and is superposed on a larger more gradual temperature trend, 
making FLI detection of the event challenging. 

D.ll Questions 

The FLI also can produce estimates of the temperature and humidity profiles by careful 
inversion of the brightness temperature relations. Does this information offer any advantages 
to identifying turbulent regions? What are the unresolved technical issues? 

The most pressing unresolved issue for evaluation of the FLI technique is the determi¬ 
nation of the false alarm rate and the detection rate for typical flight conditions. 

Do the EOF simulations include the effects of correlated noise and is this possible ? The 
large dynamic range in the eigenvalues may be sensitive to correlated noise. 
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Figure D8. Sample spatial spectra of u, v, w, and T from INDOEX 
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Figure D14. Spatial spectra of turbulence data from NOAA G4 aircraft. 
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Figure D15. Commercial transport aircraft turbulence over Iowa. 



Figure D16. Commercial transport aircraft turbulence over Tennessee. 
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Appendix E 


Airport Runway Friction Detection System 


William L. Smith 

Science System and Applications, Inc. 

George Tsoucalas 
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E.l Introduction 

Since the beginning of commercial aviation, there have been hundreds of accidents caused 
by poor runway friction induced by rain, frost, and snow. During 2007 alone, 289 people 
were killed in two accidents related to poor runway friction. On July 17, 2007, a passenger 
jet skidded off a rain-slicked runway in Sao Paulo Brazil and crashed into a cargo terminal 
and fuel station, bursting into flames. Everyone (187) on board and 13 more on the ground 
were killed in the accident, killing 200 people in the country’s worst aviation disaster. More 
recently on September 16, 2007, 89 passengers were killed when an airliner landing at Phuket 
Thailand in a torrential rain skidded off the runway, ran through a low retaining wall, split 
in two, and caught fire. Although hindsight might indicate that the pilots of these two 
aircraft made poor judgments in trying to land under such meteorological conditions, in 
most instances rain and snow do not cause the runway to be too slick for a safe landing. 
Thus, what is needed to prevent such tragic accidents in the future is a runway detection 
system, with a display in the cockpit, which can detect the runway condition just prior to 
touchdown and alert the pilot to the total length of runway needed after touchdown to come 
to a safe stop. 


E.2 Potential Airborne Sensor Approach 

The motivation for the detection system proposed here is the use of a Forward-Looking 
Interferometer (FLI) spectrometer under study for the detection of a variety of in-flight 
aviation weather hazards, including clear air turbulence. The FLI is a highly sensitive 
imaging spectrometer operating in the infrared portion of the electromagnetic spectrum. 
The ability to detect and quantify the runway friction condition with a FLI is based on the 
fact that the spectral distribution of radiance emitted by a runway within the 8-13 /zm region 
of the spectrum has a unique spectral feature depending upon whether the runway surface 
is bare concrete or asphalt, or covered with water, and/or ice, and/or snow. Moreover 
the strength of this spectral radiance emission feature, relative to that expected for a dry 
runway condition, will depend on the amount of rain, ice, or snow cover, which together 
with the known aircraft weight and landing speed, provides a quantitative measure of the 
runway coefficient of friction (i.e, COF defined as the ratio of the horizontal braking to the 
vertical loading forces acting on aircraft tires). Using an imaging interferometer to detect 
the runway spectral radiance emission enables the COF to be determined over the entire 
length of the runway surface, providing with an accurate measure of stopping distance before 
the pilot of an aircraft chooses to land. 
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E.3 Physical Basis 

The radiance observed by an imaging spectrometer takes the analytical form 

R(0) = e s T s (0)B(T s ) + [1 - e s ]r s (0) J B(T)dr(0) + J B(T)dr(0) (El) 

where all variables are wavelength (A) dependent. 9 is the angle at which the instrument 
views the ground from the aircraft, e s is the emissivity of the surface which has a different 
spectral character for different surface conditions (e.g., dry concrete, asphalt, or various 
degrees of rain, frost, snow, or ice covered surfaces), B is the Planck function, which, 
unlike the surface emissivity, is a very smooth monotonic function of wavelength. The 
transmittance of the atmosphere between the aircraft and any position along the viewed path 
is t(6), t s (6) being the transmittance of the atmosphere between the aircraft and the surface. 
The first term represents the contribution of the measured radiance coming from the surface 
and the second and third terms represent the contributions from the atmosphere. The first 
integral is the contribution of the down-welling radiation from the entire atmosphere above 
the ground that is reflected by the ground surface to the instrument, whereas the third term 
is the contribution of the upwelling radiance from the atmosphere between the instrument 
and the ground surface. For an instrument viewing the runway from an aircraft close to the 
runway surface across the 8-13 /rm “window” portion of the spectrum where atmospheric 
transmittance is close to unity, equation El simplifies to 

R(0) = e s B(T ) (E2) 

According to equation E2, the spectral character of the observed runway radiance signal 
will be dominated by the spectral character of runway surface emissivity, e s . 

Figure El shows the emissivity of various types of runway surfaces in which the runway 
surface is either dry or saturated with pure water, snow, or ice. The runway spectral radiance 
distribution, formed by a mixture of these surface conditions, will be a linear combination of 
the spectral radiance distributions forming the mixture, where the weights are the fractions 
of cover by each constituent. 

The unique spectral character of the radiance emission by each of these runway surface 
types together with the knowledge that any particular runway condition will be characterized 
by a linear combination of these spectral distributions enables one to determine the COF 
for any viewed spot under clear line of sight conditions. 

E.4 Runway Stopping Distance Retrieval Algorithm 

For a clear atmosphere between the aircraft and the runway surface, the retrieval of the 
spectral character of the surface emissivity spectrum from the observed radiance spectrum 
is straightforward. The system of linearized equations to be solved is 

R x = e x B x (T 0 ) + e x (T - T 0 ) dB ^ (E3) 

where e x is the surface emissivity, T is the true surface skin temperature, T 0 is an estimate 
of the surface skin temperature, such as the reported surface air temperature, and B x is the 
Planck function. Now assuming that 

e x = Y, a itXi,i=lA (E4) 

where a*, is the fraction of each emissivity spectrum shown in figure 1 and 

5> = l, i = 1,4 (E5) 
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Figure El. The spectral distribution of emissivity of various types of runway surfaces. 
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The combination of equations E3, E4, and E5 provide a system of N + 1 equations, where 
N is the number of spectral radiances measured, to solve for five unknowns (i.e., the four a’s 
and T). Since an interferometer spectrometer provides hundreds of spectral radiance obser¬ 
vations within the 8-13 pm region of the spectrum, the system is heavily over determined 
providing stable estimates of the unknown parameters. 

The oq’s provide an estimate of the state of the runway surface from which the COF 
can be estimated. Given the spatial distribution of COF throughout the entire length of 
the runway and given the aircraft weight and landing speed characteristics, the length of 
runway needed for stopping can be determined, dependent on touchdown point. 


E.5 Low Visibility Atmospheric Conditions 

For a low visibility atmospheric condition, which often co-exists with deteriorating runway 
surface conditions, it has been shown that there is sufficient information in the measured 
radiance spectrum to account for the attenuation of the surface emitted radiance spectrum 
for a simple constant surface emissivity condition (i.e., emissivity equal to unity). Figure 
E2 shows an example of how a poor visibility (-<1000 ft) (-<305 m) image of a runway 
obscured by haze and fog can be cleared using a Principal Component analysis surface 
radiance retrieval algorithm. The atmospheric attenuation enhanced image of a runway 
with a rectangular object in the middle using radiance observations from a Forward Looking 
Interferometer (FLI) were shown in Section ??. 



Enhanced, 8v = 0.6 cm 1 


Enhanced, 8v = 2.5 cm -1 


Figure E2. Runway object image enhancement. 

For this simulation, the instrument was assumed to be 1000 ft (305 m) above a spectrally 
“black” runway surface in which the runway temperature was one degree warmer than the 
surrounding grass and the object in the middle of the runway surface was 2 degrees colder 
than the runway. The simulation shows that the spectral character if the atmospheric 
attenuated runway emission can be used account for the optical depth of the atmosphere 
between the runway surface and the aircraft mounted instrument. However, it remains 



















to be determined whether or not there is sufficient information in a measured atmospheric 
attenuated radiance spectrum emitted from a runway with a more complex surface emissivity 
condition, such as those shown in figure El. However, assuming that the spectrum of aerosol 
and fog optical depth possesses a sufficiently different spectral and spatial character than 
the runway surface emissivity, and/or surrounding grass emissivity, it is possible that these 
two contributing components to the observed radiance spectrum can be separated given a 
sufficient observation of the three dimensional (x,y, A) radiance distribution as would be 
obtained with an imaging forward looking interferometer spectrometer. 

E.6 Conclusion 

It appears to be theoretically and technically feasible that a relatively low spectral resolu¬ 
tion (A/<5A «400) imaging FLI spectrometer, operating from an aircraft, could be used to 
detect the runway condition from which the COF and braking distance can be calculated. 
Assuming that this information could be obtained with sufficient lead-time to alert the pi¬ 
lot of a hazardous landing condition, slick runway airline disasters could be avoided. It is 
recommended that both theoretical and observational studies be performed to demonstrate 
the ability to determine hazardous runway conditions using an appropriate forward-looking 
interferometer spectrometer. 
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